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The closed, k=l, FRW model coupled to a massless scalar field is investigated in 
the framework of loop quantum cosmology using analytical and numerical methods. 
As in the k = case, the scalar field can be again used as emergent time to construct 
the physical Hilbert space and introduce Dirac observables. The resulting framework 
is then used to address a major challenge of quantum cosmology: resolving the big- 
bang singularity while retaining agreement with general relativity at large scales. 
It is shown that the framework fulfills this task. In particular, for states which 
are semi-classical at some late time, the big-bang is replaced by a quantum bounce 
and a recollapse occurs at the value of the scale factor predicted by classical general 
relativity. Thus, the 'difficulties' pointed out by Green and Unruh in the k=l case do 
not arise in a more systematic treatment. As in k = models, quantum dynamics 
is deterministic across the deep Planck regime. However, because it also retains 
the classical recollapse, in contrast to the k = case one is now led to a cyclic 
model. Finally, we clarify some issues raised by Laguna's recent work addressed to 
computational physicists. 

PACS numbers: 04.60.Kz,04.60Pp,98.80Qc,03.65.Sq 



I. INTRODUCTION 



The spatially flat, isotropic, k=0 model was recently investigated in detail in the setting 
of loop quantum cosmology (LQC) [H, Q, S]- That investigation introduced a conceptual 
framework and analytical and numerical tools to construct the physical sector of the quan- 
tum theory. These methods enabled one to systematically explore the effects of quantum 
geometry both on the gravitational and matter sectors and to extend previous results in 
LQC. The purpose of this paper is to use the 'improved dynamics' of [3] to carry out a 
similar analysis for the closed, k—1 model, coupled to a massless scalar field. In our pre- 
sentation, we will skip those constructions, proofs and arguments which are direct analogs 
of the ones used in the k=0 case in [H, 0] • The focus will be on the differences from the 
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existing treatments of the k=l model and also our earlier analysis of the k=0 model. 

Although current observations favor spatially flat models, the k=l closed model is of con- 
siderable conceptual and technical interest. On the conceptual side, an outstanding problem 
in quantum cosmology — and indeed in full quantum gravity — is whether one can construct 
a framework that cures the short- distance difficulties faced by the classical theory near singu- 
larities, while maintaining an agreement with it at large scales. By their very construction, 
perturbative and effective descriptions have no problem with the second requirement. How- 
ever, physically their implications can not be trusted at the Planck scale and mathematically 
they generally fail to provide a deterministic evolution across the putative singularity. In 
loop quantum gravity (LQG) the situation is just the opposite. Quantum geometry gives 
rise to new discrete structures at the Planck scale that modify the classical theory in such a 
way that, at least in simple models, space-like singularities of general relativity are resolved. 
However, since the emphasis is on background independence and non-perturbative methods, 
a priori it is not clear whether the theory also has a rich semi-classical sector. Do the novel 
dynamical corrections unleashed by the underlying quantum geometry naturally fade away 
at macroscopic distances or do they have unforeseen implications that prevent the theory 
from reproducing general relativity at large scales? While there is recent progress which 
indicates that LQG does admit a viable semi-classical sector near Minkowski space-time 
further evidence is needed in other contexts such as cosmological settings. 

The 'improved dynamics' of 0] successfully addressed this issue in the k=0 case. However, 
because classical recollapse, the k—1 model provides a more stringent tests. In particular 
using numerical evolutions of equations used in the early LQC works [5], Green and Unruh 
0] had concluded that there may be a key difficulty: the LQC evolution appeared not to 
reproduce the recollapse predicted by general relativity. Since curvatures at the epoch of 
the classical recollapse are very small, this feature appeared to indicate that LQC would 
deviate from general relativity in perfectly ordinary situations. Thus, although the LQC 
equations used in the Green-Unruh analysis did cure the ultraviolet difficulties by resolving 
the singularity in the sense of Q, they appeared not to have a viable infrared behavior. Can 
this difficulty be resolved? Or, does the situation in the k—1 model indicate that LQG may 
not admit a good semi-classical sector in this cosmological setting? We will employ methods 
developed in the k=0 case [3] to probe this issue in the physical sector of the quantum 
theory. A systematic analysis will show that not only is the big-bang singularity resolved 
but the quantum evolution in fact faithfully mirrors the predictions of general relativity, 
including the re-collapse, when the curvature is small compared to the Planck scale. 

A second conceptual issue is whether the general features of the Planck scale physics found 
in the k=0 case in [3J are robust. For example, there the Friedmann equation (a/ a) 2 = 
811G p/3 is replaced in LQC by the quantum corrected equation (a/a) 2 = (8irG p/3)(l — 
p/ Petit) where p crit ~ 0.82p P1 . The correction comes with a negative sign, making it possible 
for d to vanish — triggering a bounce — when the matter density reaches a critical value, 
Pcrit- One then has a pre-big-bang branch joined on to the current post-big-bang branch 
by a deterministic evolution, both behaving classically when the density is low compared to 
Pcrit- Does this feature survive in the k=l case or does space-time simply become fuzzy near 
and to the 'past' of the big-bang? Is there only one cycle in which the universe resembles 
our own? Is the value of the critical density p max at the bounce point essentially the same 
as in the k=0 model or does it depend on the spatial topology? Our analysis will show that 
the big-bang and the big-crunch are replaced by a quantum bounce leading, in a precise 
sense, to a cyclic quantum universe. Furthermore, the value of p max is robust so long as the 
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maximum radius of the universe attains a macroscopic size. 

On the technical side, this model also provides a number of challenges. In LQG, the 
configuration variable is a connection A l a , related to the spin connection T l a (determined by 
the spatial triad) and the extrinsic curvature K l a via A % a = T l a + ^K % a where 7 is the Barbero- 
Immirzi parameter. However, because of certain technical difficulties with the methods used, 
so far the k=l case has been treated in the literature by regarding the extrinsic curvature 
K l a as a connection and constructing holonomies from it [5[. Because of gauge fixing this 
is a mathematically viable strategy. However as emphasized in, e.g. |6|], to make closer 
contact with the full theory, it is desirable to construct the theory using connections A\. So 
an important question arises: Can the improved dynamics of |3| overcome the difficulties 
faced in the earlier treatments, allowing us to formulate the theory in terms of A l a l We 
will find that the answer is in the affirmative. A second technical challenge arises in the 
definition of the operator representing the field strength F l ah that features in the expression 
of the Hamiltonian constraint. In LQC, this operator is constructed from holonomies 0. In 
the k=0 model, one can use the integral curves of the (right and left-invariant) fiducial triad 
°e" to build the holonomy loops. In the k—1 case, °e" (is only left invariant and) satisfies 
the commutation relations of so(3). Hence none of the three pairs of triad vectors is surface 
forming. If one uses a general loop which is not tangential to these triads, holonomies (fail to 
be almost periodic functions of connections and) do not lead to well-defined operators in the 
quantum theory. Thus, finding suitable loops poses an interesting technical challenge. We 
resolve this issue. 1 Finally, to numerically solve the quantum Hamiltonian constraint, one has 
to find normalizable eigenfunctions of a certain difference operator 0. Recall that in the case 
of, e.g., a simple harmonic oscillator, the differential operator representing the Hamiltonian 
admits an eigenfunction with any real number as its eigenvalue. Normalizable eigenfunctions 
exist only for discrete eigenvalues and thus constitute a set of zero measure in the space of 
all eigenfunctions. In the k=l case, the situation is similar with the operator 0. However, 
now the normalizable eigenfunctions and their eigenvalues are not known analytically and, 
since they constitute a set of zero measure, one has to exercise considerable care in singling 
them out numerically. 

Insights gained from the resolution of these conceptual and technical issues are likely to 
be important as one considers more and more general situations and develop new tools to 
tackle the Hamiltonian constraint of full LQG. 

The paper is organized as follows. In section [TT] we summarize the basics of the classical 
and quantum theories for our model, arriving at the form of the Hamiltonian constraint op- 
erator. In section HTT1 we sketch the Wheeler-DeWitt (WDW) theory and in II VI we introduce 
the physical sector of the theory in LQC. Numerical simulations are discussed in section 
|Vl In section ED we summarize effective equations of the semi-classical sector, list the main 
numerical results and compare them with predictions of the effective theory. We also clarify 



some issues that have been raised by a numerical analysis by Laguna and remove confu- 
sion caused by an unfortunate use of terminology there. Section IVHl places these results in 
a broader context and discusses their relation with that by Green and Unruh ||. Appendix 
IA1 summarizes some background material on the fiducial geometrical structures used in k=l 
models with the so-called 'Bianchi IX' symmetries. 



1 This issue was resolved independently by the Warsaw group (lo| . Their method of evaluating the holonomy 
is more intrinsic, elegant and insightful. 
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II. THE BASIC FRAMEWORK 

This section is divided into two parts. In the first we summarize the classical theory and 
in the second we present the basics of the quantum theory. 



A. Classical Theory 

Space-time manifolds under consideration are of the form M x K, where M has the 
topology of a 3-sphere, S 3 . As explained in appendix E] one can identify M with the 
symmetry group SU(2) (which ensures spatial homogeneity and isotropy) and endow it with 
a fixed fiducial basis of 1-forms °uo l a and vectors °e\. The resulting fiducial metric is 

°Qab '■= 0uJ a 0uJ l kij' the Cartan-Killmg metric on su(2). (2.1) 

°q a b turns out to be the metric of the round 3-sphere with radius a Q = 2 (rather than 
a Q = 1). The volume of M w.r.t. this fiducial metric °q a b is V Q = 2ir 2 a? = 167r 2 and the 
scalar curvature is °R = Q/a 2 = 3/2. We shall set i Q := Vl^ . (For details, see Appendix [AT) . 

To facilitate comparison with the spatially flat, k=0 case and to directly use results 
of [l|, 0, S] obtained in that case, we will set our conventions in a parallel fashion. The 
dynamical, isotropic homogeneous connections and triads will be parameterized by c and p 
respectively: 

A\ = cl- Xo ^ Et = pC^q°e1. (2.2) 

As in the k=0 case c is dimensionless while p has dimensions of area and the density weight 
of Ef is absorbed in the determinant of the fiducial metric. At the point (c, p) of the phase 
space, the physical 3-metric q a b and the extrinsic curvature K ab are given by: 

Qab = \p\ C\ab and 1 K ah = (c- f ) \p\\ £- 2o q ab (2.3) 

3 

The corresponding physical volume of M is \p\ 2 - The scale factor a associated with a 
physical metric q ab is generally expressed via q ab = a 2 °o Lab where °q ab is the unit 3-sphere 
metric. Then, the scale factor is related to p via \p\ = a 2 l 2 /A. As usual p can take both 
positive and negative values, the change in sign corresponds to a flip in the orientation of 
the triads °e" which leaves the physical metric q ab invariant. 

Expressions of the total action, the symplectic structure, and the Hamiltonian constraint 
require an integration over M. In the k=0 case the underlying manifold is non-compact 
(1R 3 ), whence the naive integrals would simply diverge. Therefore, irrespective of which 
quantization scheme one uses, one has to fix a fiducial cell with finite volume, say V (w.r.t. 
the fiducial flat metric), and restrict all integrations to this cell. In the k—1 case, M is 
compact and the introduction of a cell is unnecessary. Nonetheless, in many of the key 
equations our V plays the role of the volume V of the cell in the k=0 case. In both cases, 
the fundamental Poisson bracket is: 

{c, P} = (2.4) 



where 7 is the Barbero-Immirzi parameter. Finally, using the fact that the Cartan orthonor- 
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mal triad u l a on a 3-sphere of radius a Q satisfies the identity 

du k + — eifu 1 A u j = , 



(2.5) 



it is straightforward to calculate the field strength F k b of the connection A l a on M 

. , 2 



1 ab c o 



c 2 — C\ 



_ k o j q k 
€ ij U a U b ■ 



(2.6) 



Our fiducial co-triad °uj % a corresponds to a Q = 2. But we have refrained from using this 
numerical value in (12.61) to clarify the relation with the k=0 case. To pass to this case, one 
has to set V a equal to the volume V D of the fiducial cell there but take the limit as a Q goes 
to infinity. In this limit, the fiducial co-triad °uj l a goes to the (flat) fiducial co-triad used in 

and the field strength F* b goes over to its value V Q 2 ^c 2 €ij k 



the k=0 case (see (12.51) ) 
in the k=0 case [H, 0, H 

As in the k=0 models, we have completely gauge fixed the Gauss and the diffeomorphism 
constraints by fixing °io % a and °e" and using the parametrization (12.21) of the phase space 
variables A l a , E" 1 . So, we are left only with the Hamiltonian constraint. The gravitational 
part of this constraint is given by [l2j, [l5| : 



a 



grav 



d 3 x 



M 



E ai E bj 



rpk 

fab - 



2(1 + 7 2 ) 



E ai E bj 



(2.7) 



where e := y^jdetEj and K % a = K a ho ujl is the extrinsic curvature, and where, as is usual in 
mini-superspace analysis, we have set the lapse equal to 1. Using the fact that A\ is related 
to the spin-connection T l a (of the physical triad ef) and the extrinsic curvature K l a through 
A l a = T l a + / -fK] 1 , it is convenient to express the second term in the integrand of (12.71) in terms 



of the curvature F* b of A\ 



and the curvature of T l 



E ai E bi K i aR ^ = 



2r 



i E bj ^ 



k 
ab 



^ab) 



(2.8) 



In the k=0 case, the spatial curvature vanishes and the extrinsic curvature term in (12. 7p 
reduces to a multiple of the first term involving F k b . In full general relativity, by contrast, 
while Q l ab is determined by the momenta Ef , its expression is rather complicated and this 
strategy of handling the extrinsic curvature terms, by itself, does not simplify matters. The 
situation in general homogeneous models [E| as well as black hole interiors [l3[ is in between 
the two: Although Q k b is non-zero, its expression is simple and manageable. In the k=l 
model now under consideration one has: 



Qk 



1 



°l !* °, J 



(2.9) 
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Therefore, the gravitational part of the constraint reduces to: 



a 



grav 



1 

7 2 JM 



* 3 



ab 



{ 4 ) e "b U c 



r 



{ 2 aj 4 o? n 



(2-10) 
(2.11) 



where, in the last step we have used the expression ( 12.61) of F a f,. In the k—1 case now under 
consideration, a Q = 2. However, as indicated above, we did not substitute this numerical 
value because results for the k=0 case [l|, 0, [3| can be recovered by setting a Q = oo (and 
£ = V )• However, since a Q always occurs in the combination £ /a in the Hamiltonian 
constraint, we can set a D = 2 throughout and recover the k=0 results simply by setting 
£ = 0. In what follows, we will adopt this strategy. 

Remark: In the above construction, we began with a fiducial triad °ef and a co-triad °uj l a 
adapted to a 3-sphere of radius a Q =2. Therefore, our construction may appear to be tied to 
that choice. Had we used a 3-sphere of radius a Q = 2A, the fiducial triad and the co-triad 
we have rescaled via °e" — > A -10 ^ and O uj^\°oj l a . It is easy to check that the variables c, p 
parameterizing the physical fields (A l a , Ef) are left unchanged. Hence the entire framework 
is invariant under this rescaling 'gauge' freedom. 



B. Quantum kinematics and the Hamiltonian constraint 

To pass to the quantum theory, following Dirac one first constructs a kinematical de- 
scription. As in the k=0 case the kinematical Hilbert space TCf^ is the space 
L 2 (MBohr, ^Bohr) of square integrable functions on the Bohr compactification of the real 
line. To specify states concretely, it is convenient to work with the representation in which 
the operator p is diagonal. Eigenstates of p are labeled by a real number /i and satisfy the 
orthonormality relation: 

(Mi|M2>=<W 2 - ( 2 - 12 ) 
Since the right side is the Kronecker delta rather than the Dirac delta distribution, a typical 
state in can be expressed as a countable sum; = J2 n ^ll^n) where are complex 
coefficients and the inner product is given by 

(* 1 |* 2 ) = 53^ ) c4 n) • (2.13) 

n 

The fundamental operators are p and expiA(c/2): 

P\fi) = — and exp — = |/i + A) (2.14) 

where A is any real number. From the discussion of the classical theory of section III Al it 
follows that the physical volume operator of M is given by: V = \p\ 3 ^ 2 . 

Of special interest to us are holonomies of the connection A l a along the integral curves 
of our fiducial triads °ef. The holonomy along the segment of (directed) length X£ , 
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tangential to e\ is given by: 



.2 



hP = cos — 1 + 2 sin — r k . (2.15) 



The corresponding holonomy operator has the action: 

h[ X) I//) = \ (I// + A) + |/i - A)) I + \ (I// + A) - |/i - A» r k . (2.16) 
However, just as there is no operator corresponding to the connection itself in full LQG 

/grav 
L kin 



121 . 1151 ] , there is no operator c on |8j 



To describe quantum dynamics, we have to first introduce a well-defined operator on 
representing the Hamiltonian constraint C grav . Since there is no operator corresponding to 
c itself, as in the k=0 case we will use the integral expression H 2 . 101) of the constraint. For 
the passage to quantum theory, one has to first express this classical constraint in terms of 
the elementary variables p and and then replace them with operators p and ^As in 



the full theory [14], [15J , the term involving triads becomes the following operator 



e ijk er^W* = £ jTj^ ^ Tr [A^" 1 , V] n) (2.17) 

k 

where V = \p\ 3 ^ 2 is the volume operator. 

To define F ab k in f)2.10p . as in [l|, 0, H, [1] we use the standard relation between the 
holonomies and field strengths. Because of homogeneity, the components °e" °e b j F a b k are 
constant on M. They can be evaluated by considering a square loop starting and ending 
at any point x, with tangent vectors °ef and °e b j at x, and then taking the limit as the area 
enclosed by the loop shrinks to zero. In the quantum theory, since there is no operator 
corresponding to the connection c, the limiting operator does not exist. As discussed in 
detail in [j], @, H, 0], this is a manifestation of the quantum nature of geometry, i.e., a 
reflection of the fact that the area operator has purely discrete spectrum. As in |3( our 
strategy is to shrink the loop only until its physical area equals the 'area gap' A, i.e., the 
minimal non-zero eigenvalue of the area operator. 

Now, in the k=0 case, the edges of can be taken to be the integral curves of the 
triad vector fields °e" and °e^. In the k—1 model, however, the °e" satisfy the commutation 
relations of su(2); they do not commute. Therefore their integral curves can not provide the 
desired closed loop Dy. In the existing literature, a closed loop is formed by simply adding 
a fifth edge [1]. However, this strategy is not viable: Since the five edges do not span an 
unambiguous 2-surface, the notion of the area enclosed by the loop has no obvious meaning. 
However, note that while °e" are the 'left-invariant' vector fields, M also admits three 'right 
invariant' vector fields £ 4 a (see Appendix |A]) . These are the symmetry fields: they also satisfy 
the commutation relations of su(2), act simply and transitively on M and their action leaves 
each of our fiducial triads °e° and co-triads to^ invariant. Since they commute with °e", one 



2 Here I is the unit 2x2 matrix and is a basis in the Lie algebra su(2) satisfying t, Tj = \tijhT k — 
Thus, 2iTk — (Tfe, where <Ji are the Pauli matrices. The directed length is positive if the line segment is 
oriented along e£ and negative if it is oriented in the opposite direction. 
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can form the desired closed loops by first moving from x along the integral curve of say 
%i then £k then along — °ej and then along — (where is chosen to coincide with °ej at x). 
An explicit realization of this procedure is presented in Appendix |A] and a more geometric 
construction appears in The final field strength operator does not depend on whether 
the first segment is chosen to be left invariant or right invariant and is given by 3 : 



ab 



lim 



Ar Dij-*A 



(sin 2 A(c - 4/2) 



sin 



(sin 2 fx(c - £ /2) - sin 2 (/2£ /2) ) e, 



Here, as discussed in detail in [3|, // is a specific function of p 



f \p\ = A = (2V37T 7 ) 



pi > 



(2.18) 



(2.19) 



and for notational simplicity we have dropped hats on operators which are trigonometric 
functions of c. The fact that ft is a function of p rather than a constant arises from the 
requirement that the physical area of be set equal to A. As explained in [3], this strategy 
mimics the full theory in a well-defined sense and the resulting 'improved dynamics' is free 
of the drawbacks of older Hamiltonian constraint of LQC As in [H, Q, 0, @] , the viewpoint is 
that at the fundamental level — i.e. at the Planck scale — the field strength operator is non- 
local and the usual local classical expression arises only on coarse graining in semi-classical 
states. 

To obtain the explicit action of F ab k on , one has to face two complications. The 
first arises already in the k=0 case. The operator F a f ) k for k=0 can be recovered by setting 
£ Q = 0. It depends on the connection only through sin(/Ic). Since Jx is itself a function of p, 
the action of sin llc on 7^^ v is rather subtle. As discussed in detail in [3] , it is simplest to 
express it by going to a basis \v) which is better adapted to the volume operator V: 



V\v) 



7r l\^ \ V \ b3 l \ 

— tpi \V) 



6 



K 



(2.20) 



where the dimensionless label v — the eigenvalue of V apart from a constant — is related to 
the dimensionless label \l — the eigenvalue of p apart from an overall constant — via 



2v/2 
3\/3\73 ' 



v = Ksgn(Li) |/i| 2 , where K 

In this basis, expiflc are simply the translation operators: 

e iJ f = + 1), 

so that 

sin(/Zc) = — + 2) - - 2) ). 

— t 



(2.21) 

(2.22) 
(2.23) 



3 Here and in what follows, in light of results of 
tation. For a more detailed discussion, see 0, |i|. 



16|, [17[, we have used the fundamental, j = 1/2 represen- 
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In the k=l case we have an added complication: (12.181) contains sin fi{c — i /2) rather 



than sin(/xc). This difference can be handled as follows: On 7i 



grav 
kin 



sin(/ic - ^) V(v) = e u ° f sin/xc e~ li ° f y{v), where / = 1-^1*. (2.24) 

2 4 K 

Note that since /(f) is continuous everywhere (including the point v — 0), the operator 
expi£ f is well-defined and unitary on all of . 

We now have operators corresponding to each term in the integrand of the gravitational 
part fl2TT0D of the Hamiltonian constraint. Using f[2TTTj) . fl2~T8|) . f[2T23l and fl2T24l in fl2TT0|) 
we obtain: 



C'grav^(w) = sin /ic A sin /2c e~ i£o/ ^ 

~2 4 

where, as in the k=0 analysis (3| we have set 



[ sin 2 ^ - — 2 ] A ^(f ) (2.25) 

9|A 2 d3 



^^) = -^Vy^H |l«-l|-k + l|| *(«)■ (2.26) 

To summarize, there are two main subtleties in the passage from spatially flat, k=0 models 
to the closed, k=l ones. First, the loop around which holonomy is computed has to be 
constructed using both right and left invariant vector fields. Second, now the connection 
dependence is in the operator sin/2(c — £ D /2) rather than sin flc and one has to define its 
action on 7if^ v using the unitary operators exp i£ a f. When this is done, the gravitational 
part of the Hamiltonian constraint is symmetric and positive, given by (12.251) . 

Finally, to write the complete constraint operator we also need the matter part of the 
constraint. For the massless scalar field, in the classical theory it is given by: 

C matt = 87rG \p\~lpl (2.27) 

As usual, the non-trivial part in the passage to quantum theory is the function |p|~ 3 / 2 . How- 
ever, as with the co-triad operator ( 12.171) . we can use the method introduced by Thiemann 
in the full theory 

HI 

Il5| . This issue is discussed in detail in [Ijj]. The final result is: 



3/2 



\p\-2^( v ) = ( — — ) B{v) V(v) (2.28) 



where 

B(v) =(-) A' H 



3 x3 



+ l|l/3_| l; _ 1 |l/3 3 . 

It is self-adjoint on 7ij^ v an< ^ diagonal in the eigenstates of the volume operator. 
Collecting these results we can express the total constraint 



(2.29) 



C ^{v) = ((7 grav + C^) («) = , (2.30) 
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as follows: 



1 

3 



— [Biv)}- 1 3K(sin 2 (^) Jj2f °- 



2 



St 2 













1 \v — 1| — \v + 1| 


1 1 . 





2 



4 



(2.31) 



Here, o is the operator that appears in the k=0 quantum constraint in place of 0]: 

9 *(«, 0) = -[Biv)]- 1 (C + (v) *(v + 4, 0) + C» 0) + C"» *(v - 4, 0)) (2.32) 
where the coefficients C ± (f) and C°(V) are given by: 

3ttKG 



\v + 2\ \\v+l\ 



v + 3 



C~(v) = C + {v-A) 
C°(v) = -C + (v)-C-(v) 



(2.33) 



Thus, the k—1 quantum constraint has the same form as in the k=0 case. As one would 
expect from the classical expression (12.101) . the difference — O is diagonal in the v- 
representation and vanishes when we set i a = 0. 

In the remainder of the paper, we will work with the Hamiltonian constraint (12.311) . As in 
the k=0 case @, 0], the form of this constraint is similar to that of a massless Klein-Gordon 
field in a static space-time, but now with an additional static potential. is the analog of the 
static time coordinate and the difference operator 0, of the spatial Laplace-type operator 
plus the static potential. Hence, the scalar field can again be used as 'emergent time' in 
the quantum theory. We will examine the operator in some detail in sections [IV] and |Vj 
Finally, in the above construction we made a factor ordering choice the motivations behind 
which are the same as those discussed in [^]. This choice will facilitate comparison between 
the LQC results in the k=l and k=0 cases and yield the WDW equation with its 'natural' 
factor ordering in the 'continuum limit'. 



Remark: Velhinho [18( has pointed out that in quantum kinematics it would suffice to 



consider the algebra generated by p and just two almost periodic functions of the connection, 
e i^ic anc j e ifi 2 c^ w h ere n x / n 2 is irrational, because these functions already separate points of 
the phase space. In the older, '/i G -e volution' 0, [1], this strategy would be natural; one 
could set Hi — 1 and /i 2 = yU = 3\/3/2. However, as discussed in 0], this evolution is 
not viable physically. The 'improved' /i-evolution used in [3| and in this paper is free of 
those drawbacks. But since p, is now not a constant but a function (12.191) of the Velhinho 
kinematics will not support the 'improved quantum dynamics.' 



III. WHEELER DEWITT THEORY 



In this section we will briefly discuss the WDW limit of LQC in which effects specific to 
quantum geometry in the difference equation fl 2 . 3 1 [) are ignored by letting the area gap go to 
zero. This discussion will bring out the key role played by quantum geometry near the big 
bang and the big crunch singularities. The WDW theory has its roots in geometrodynamics 
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which is insensitive to the choice of the triad orientation — i.e., to the sign of v. Therefore, 
as in @, 0] , we will restrict ourselves to wave functions }&_(v) which are symmetric under the 
orientation reversal operator II, 

n V(v, 0) = 0) . (3.1) 

where (and in what follows) we have denoted the WDW analogs of the LQC quantities with 
an underbar. As in [3j we will find that the scalar field can again be used as emergent 
time and, in the resulting physical sector of the theory, the big bang and the big crunch 
singularities persist in the WDW limit. Details of motivation and background material as 
well as underlying assumptions and technical steps involved in the WDW limit can be found 




A. The WDW constraint and its general solution 

To pass to the WDW theory, one sends the area gap to zero and restricts oneself to smooth 
wave functions }&,(v,<j)). As explained in 0, 0], this corresponds to taking the continuum 
limit of the difference equation (I2.3ip . The non-trivial part of the task lies in the limit of O 
since the remainder is diagonal in v. This has already been completed in (3|. The required 
WDW limit is then given by 4 

9l%{vA) = -my A) 

= -QMvA)-— -£h«|»£M), (3.2) 

SK 3 

where, as in [§], 

© *K 0) := -12ttG (vd v ) 2 0) . (3.3) 

As explained in [3], this factor ordering is 'covariant' from the geometrodynamical perspec- 
tive and coincides with the one used in the older WDW literature (see e.g. [191]). 

The operator is self-adjoint and positive definite on the Hilbert space L 2 (M, B(i> )dv ) = 
L 2 (M, (K/\v\)dv), where, as in 0], to facilitate comparison with LQC results we have denoted 
the WDW limit K/\v\ of B(v) by B(i>). The general solution of (13. 2p can be readily expressed 
in terms of the spectral family of 0. Let us begin by considering all eigenf unctions of the 
differential operator 0: 

- 12ttG {vd v f%{v) + ^ ° 4 7 \v\^ u (v) = lu 2 Mv). (3.4) 
In their most general form, they can be expressed as linear combinations of modified Bessel 



4 The appearance of the Barbero-Immirzi parameter 7 in the WDW limit is an artifact of our conventions, 

i.e., definition of v. If instead we use the eigenvalue v of the volume operator, V\v) — vi^ \v), so that 
3 

v = (8717/6)2 (v/K), the WDW equation (|3.2[) would be manifestly independent of 7: dv$(v,(f)) = 
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functions X. /C 19l 



lM«) = «/C, fc (-^Hl) + /3 l ik (^\ v \l) , (3.5) 
V4K3 / \4i^3 / 

1 

where k := (3/lQnG) 2 u and a, (3 are constants. (We have used calligraphic letters to 
denote the modified Bessel functions to avoid confusion with the constant K of (12.211) ). As 
x — > oo, the function X ik {x) grows exponentially whereas /C^(x) decays exponentially. 

These properties of eigenf unctions imply that the spectrum of the operator is con- 
tinuous and that Iik{x) cannot feature in its spectral decomposition. Therefore to obtain 
this decomposition, we will set (3(k) =0 and choose constants a(k) such that the resulting 
eigenfunctions e k (v ) := a(k)JCik(v) are orthonormal, 

(e fc |ej.,) = 8(k,k') . (3.6) 

The Dirac distribution appears on the right hand side because the spectrum of G is continu- 
ous. This may appear surprising because e k (v) decay for large v. However, e k (v) fail to have 
finite norm in L 2 (IR, (K/\v\)dv) because they have the following oscillatory form for small v 

r f t-ol , | A |«|<* ( v / u . , , ( A ( _s 

rLik \ 2"M 3 > a \ UJ ) cos = in p| + cr(a;j , (3.7) 

\AKz J \Vl2nG J 

where cr(u) is a constant (but to- dependent) phase shift. 

In the spatially open, k=0 models [3j, the spectrum of 6 is 2-fold degenerate, reflecting 
the fact that in the classical theory there are two sets of distinct solutions, one perpetually 
expanding and the other perpetually contracting. In the present case, this degeneracy is 
broken because T ik {v ) diverge at large v. Its classical counterpart is the fact that, because 
each solution has both contracting and expanding epochs, we no longer have two distinct 
sets of universes. 

Eigenfunctions e k {v) provide the standard spectral decomposition on L 2 (R,B(v)dv): 

POO 

= / dk^(k)e k (v). (3.8) 
Jo 

Hence, the general solution to (13.21) with smooth initial data consisting of rapidly decreasing 
functions can be written as 

/■OO 

£(v,0)= / dki> + {k)e k {v)e 1 ^ + ^>-{k)e k {v)e- iul4 ' (3.9) 
Jo 

for some suitably regular function \&±(fc). Following the terminology generally used in the 
Klein-Gordon theory, the solution will be said to be of positive (resp. negative) frequency 
if ty_(k) (resp. ty+(k)) vanishes. 

As usual, the positive and negative frequency solutions satisfy first order 'evolution' 
equations, obtained by taking a square- root of the constraint (13.21) : 

=F % d+Wiy, 0) = ^6*0, 0) . (3.10) 
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If f(v) is the initial data for these equations at 'time = O ', the solutions are given by: 

%±(v, <p) = e^V® f(v) . (3.11) 



B. Physical sector of the WDW theory 

Solutions (13.91) to the WDW equation are not normalizable in 7i™in W (because zero is in 
the continuous part of the spectrum of the WDW operator). Therefore, one has to use one 



of the standard methods |20|, [21[ to construct the physical Hilbert space H^y- Since the 
procedure is completely analogous to that used in [2, |3|, we will simply summarize the final 
results. 

Hphy consists of positive frequency solutions ^_(v, 0) to (13.21) which are symmetric under 
the reversal of the orientation of the triad, i.e. satisfy ^_(v, (p) = }&_(—v, (p), and which have 
a finite norm w.r.t. the inner product: 

(^il^ 2 ) P hy= / dvB(v)Mv,<P)%(v,<P) (3.12) 

where for notational simplicity here (and in what follows) we have dropped the subscript 
+ denoting positive frequency. On this space, a useful complete set of Dirac observables is 
provided by the momentum p^ of the scalar field, 

P**M):=-tt^^ (3-13) 
and the operator |i>L corresponding to volume at the emergent time <fi = <p , 

1% <P) = ^ M <Po). (3.14) 

Using the physical Hilbert space and this complete set of Dirac observables we can now 
introduce semi-classical states and study their evolution. Let us fix an 'instant of time' 
4> = (fi and construct a semi-classical state which is peaked at p^ = p^ and |v|^ = v*. Since 
we would like the peak to be at a point that represents a large classical universe, we are led 
to choose v* ^> 1 and p^^> H (in the natural classical units c—G—1). The second condition 
is necessary to ensure that the universe expands out to a size much larger than the Planck 
scale. At 'time' (ft = <f) , consider the state 

f°° ~ ~ (k-k*) 2 

V(v,<f> ) = dk *(jfe) e k (v)e luJ ^°- r \ where *(jfe) = e — . (3.15) 
Jo 

Here 



k* = Pi and cf>* = <\> Q + ^— cosh" 1 [(^^) 1/2 (*T 2/3 ] ■ (3-16) 

In the spatially flat case, the eigenfunctions ek(v) were just plane waves [3] and one could 
evaluate this integral analytically. The modified Bessel functions are much more compli- 
cated. Therefore, in the closed model we have to use numerical methods. They show that 



14 
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FIG. 1: a) The absolute value of the WDW wave function. The initial state at 'time (fr was 
constructed using f)3. 15j) and is sharply peaked at p^, and v*\ ( j >o . It remains sharply peaked on 
the classical trajectory with p^ = p^, passing through v* at <p = <^o- For clarity of visualization, 
only the values of |^| greater than 10 -4 are shown. Being a physical state, is symmetric under 
v — > —v. In this simulation, the parameters were: p^ = 5000, and Ap^/p^ = 0.02. b) The 
expectation values (and dispersions) of | | ^ are plotted for the WDW wave function and compared 
with the classical trajectory. The WDW wave function follows the classical trajectory into the 
big-bang and big-crunch singularities. 



this state is sharply peaked at values v*, of our Dirac observables (see Fig. [T>). 

We can now ask for the evolution of this state. Does it remain peaked at the classical 
trajectory defined by = p*, and passing through v = v * at <fi = O ? This question is easy to 
answer because (13.111) implies that the (positive frequency) solution \I/(t> , 0) to (13.21) defined 
by the initial data (I3.15P is obtained simply by replacing O by <fi in (I3.15p ! Since a, the 
measure of dispersion in (I3.15p . does not depend on 0, it follows that ^_(v , <p) continues to be 
peaked at a trajectory which is precisely the classical solution of interest. This is just what 
one would hope during the epoch in which the universe is large. However, the property holds 
also in the Planck regime and the semi-classical state simply follows the classical trajectory 
into the big-bang and big crunch singularities. (See Figs. [T^, anddb.) In this sense, the 
WDW evolution does not resolve the classical singularity. 

We will show in the next three sections that the situation is very different in LQC This 
can occur because the WDW equation is a good approximation to the discrete equation 
only for large v. Furthermore, as discussed in Q], the approximation is not uniform but 
depends on the state: in arriving at the WDW equation from LQC we had to neglect \I/ 
dependent terms of the form 0(v n ~ 3 ^^) for n > 3. For semi-classical states considered 
above, this implies that the approximation is excellent for v ^> k* but becomes inadequate 
when the peak of the wave function lies at a value of v comparable to k*. Then, the LQC 
evolution departs sharply from the WDW evolution. We will find that, rather than following 
the classical trajectory into the big bang singularity, the peak of the LQC wave function 
now exhibits a bounce. Since large values of k* are classically preferred, the value of v at 
the bounce can be quite large. However, as in the k = models [§], we will find that the 
matter density at the bounce point is comparable to the Planck density, independently of 
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the precise value of k* so long as ^> H. 

IV. LOOP QUANTUM COSMOLOGY: ANALYTICAL ISSUES 

Since generic solutions in the closed, k=l cosmologies undergo a recollapse in classical 
general relativity, the scale factor can not serve as a global time variable. This fact has been 
used as a criticism of the implicit use of the scale factor as time in the older LQC literature 
(see, e.g., 0). However, if matter sources include a massless scalar field 0, that field is 
monotonic and single valued in all classical solutions. Therefore it can be chosen as the 
'internal' time variable with respect to which the scale factor and other fields evolve. Fur- 
thermore, as discussed in section [Til f° r the model considered in this paper the Hamiltonian 
constraint (12.31j) of LQC has the same form as the wave equation in a static space-time, 
with <p playing the role of time and of the elliptic spatial operator plus a static potential. 
Therefore as in the k=0 case (f> can be regarded as emergent time also in the quantum theory. 
This choice is free of the criticisms mentioned above. 5 In this section we will construct the 
physical sector of LQC by exploiting this fact. As mentioned in section [TJ we will find that 
quantum geometry effects resolve both the big bang and the big crunch singularities, leading 
to a cyclic quantum universe. Still the scalar field will continue to be a viable emergent time 
globally. 

A. General solution to the LQC Hamiltonian constraint 

As in the spatially flat k=0 case, our kinematical Hilbert space will be: Ti-^ 1 '■= 
L 2 {R Bohr , B(v)dfiBohr) <£> L 2 {R, d(f>). Since <p is to be thought of as 'time' and v as the gen- 
uine, physical degree of freedom which evolves with respect to this 'time', we have chosen the 
standard Schrodinger representation for but the 'polymer representation' for v to correctly 
incorporate the quantum geometry effects. This is a conservative approach in that the results 
will directly reveal the manifestations of quantum geometry. Had we chosen a non-standard 
representation for the scalar field, these effects would have been mixed with those arising 
from an unusual representation of 'time evolution' and, furthermore, comparison with the 
WDW theory would have become more complicated. 

The form of the LQC Hamiltonian constraint (12.311) is the same as that of the 
WDW constraint (13. 2p . Properties of G D analyzed in [3J] immediately imply that G is again 
a positive, symmetric operator on L 2 (R Rntir , B(y)diinnhr), whence it admits a self-adjoint 
(Friedrich) extension. (For precise domains, see [lOj.) The main difference is that while the 
WDW Q is a differential operator, the LQC G is a difference operator. This gives rise to 
certain technically important distinctions. For, now the space of physical states — i.e. of 
appropriate solutions to the constraint equation — is naturally divided into sectors each of 



5 Furthermore, as emphasized in [2|, ua], while the availability of a global time variable simplifies the con- 
structions considerably and makes physical interpretation transparent, it is not essential. Using the group 
averaging procedure [2l|, for example, the physical sector of the theorycan be constructed even when a 
global intrinsic time does not exist. Furthermore, in simple examples |22| physical interpretation can be 
aided by the introduction of a suitable local notion of time which exists, e.g., if the scalar field comes with 
an inflationary potential. 
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which is preserved by the 'evolution' and by the action of our Dirac observables. Thus, there 
is super-selection. Let C\ £ \ denote the 'lattice' of points + 4n, n G Z} on the f-axis, C-\ e \ 
the 'lattice' of points {— |e| + An, n G Z} and let C £ = C\ £ \ U £-| e | where as usual Z denotes 
the set of integers. Let Wy av , WfigT and Hf av denote the subspaces of L 2 (R Bo hr, B(v)dfiBohr) 
with states whose support is restricted to lattices C\ £ \,£-\ e \ and C £ . Each of these three 
subspaces is mapped to itself by which is self-adjoint and positive definite on all three 
Hilbert spaces. However, for reasons explained in detail in [3j, our physical states will be 
invariant under the orientation reversing operator II of (13.11) . Thus, we are primarily inter- 
ested in the symmetric subspace of 7if av ; the other two Hilbert spaces will be useful only 
in the intermediate stages of our discussion. 

Our first task is to explore properties of the operator 6. Since it is self-adjoint and 
positive definite, its spectrum is non-negative. Therefore as in the WDW theory we will 
denote its eigenvalues by uj 1 . On each Hilbert space WS*T, we can solve for the eigenvalue 
equation Qi/) u (v) = J 1 ip^v ), i.e., 



+ 



C + (v) Mv + 4) + C°(v) Mv) + C-(v) Mv - 4) 
— [3X(sm (— ) - _) H - -t ol |-| 

lu 2 B(v)Mv). (4.1) 



\v + II 



Since this equation has the form of a second order recursion relation, as in the k=0 case one 
might expect a 2-fold degeneracy. However, there is an important subtlety. Let us consider 
the asymptotic regime ))>1. Then, each i]) u approaches a solution to the WDW equation 
and as we saw in section IIII Al there is only one linearly independent solution which does 
not diverge for large v. The form of the inner product on 7Y| rav now implies that, as in the 
WDW theory, the degeneracy is removed; only one of the two linearly independent solutions 
can belong to the spectral family of O. Furthermore, numerical calculations show that 
this solution, which decays exponentially for v ^> 1, will in general diverge in the other 
asymptotic limit — v 3> 1. It is only for some discrete values uo n of uj that one obtains 
solutions which do not diverge in either asymptotic limits. (This phenomenon was already 
noted in [9| in a simpler model without a scalar field but with a negative cosmological 
constant.) Since these solutions decay exponentially in both limits, they are normalizable 
in ?i| rav . Thus, in contrast to the k=0 case, on each superselected sector, the spectrum of 
O is discrete and each of its eigenvalues is non- degenerate. 6 We will denote the normalized 
eigenfunctions in TC^ V by e n {v): 

O e n (v) = u 2 n e n (v) and (e m | e n ) = 5 mi „ . (4.2) 

Finally, physical states will be built from eigenfunctions which are symmetric under 
orientation reversal. Since O commutes with the orientation reversal operator n, e n (— v) = 



See [10( for an analytical proof which does not rely on numerical results. Recall that in the WDW theory 
the spectrum is non-degenerate but continuous. This difference arises because in the WDW theory the 
eigenfunctions oscillate more and more wildly as v goes to zero and thus fail to be normalizable in the 
WDW Hilbert space L 2 (]R, (K/v)dv). Because of the discreteness of lattices C±\ e \ , nothing special happens 
'near' v = in LQC. 
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Ue n (v) is an eigenfunction of 9 with the same eigenvalue uj^ as e n (v), but belongs to W5*T 
rather than 7iy av - Therefore, 

e f) = -L( CB («) + efl (- w )), (4.3) 

also has eigenvalue cu^, but belongs to TCf 3 ". 

We can now write down the general symmetric solution to the quantum constraint (12.3ip 
with initial data in 7if rav : 

0) = K ^ + e< 8) 00 (4-4) 

where are square-summable. As in the WDW theory, if *f>~ vanishes, we will say that 
the solution is of positive frequency and if vanishes we will say it is of negative frequency. 
Thus, every solution to (12.311) admits a natural positive and negative frequency decompo- 
sition. The positive (respectively negative) frequency solutions satisfy a Schrodinger type 
first order differential equation in <p\ 

Ti-^ = Ve* ± (4.5) 

with a Hamiltonian \/0 (which is non-local in v) . Therefore the solutions with initial datum 
\I/(t>, 4> ) = f±{v) are given by: 

*±M) = e Hy/S ^ M f±(v,(f>). (4.6) 

To summarize, the overall structure is analogous to that in the spatially open, k=0 
case. A key difference is that the spectrum of is discrete on each of the three Hilbert 
spaces, TI^TZ and 7i| rav . In particular, while in the k=0 case all eigenfunctions have an 
oscillatory asymptotic behavior, now they all decay exponentially for sufficiently large \v\. 
This difference neatly encodes in the quantum theory the key qualitative difference between 
the two models in the classical theory: In the k=0 case any one classical solution is either ever 
expanding or ever contracting, while in the k—1 case each solution expands to a maximum 
volume and then recollapses. 



B. The Physical sector 

We will now summarize the mathematical structure of the physical sector of the theory. 
The construction is entirely analogous to that in the k=0 case since the spectrum of the full 
constraint operator + is still continuous (because of the <9^ part). Therefore we will 
only state the final results. 

The sector of the physical Hilbert space Hp hy labeled by e G [0, 2] consists of positive 
frequency solutions $?(v,(j)) to H4.5|) with initial data ^f(v,(j) ) in the symmetric sector of 
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^grav Eq. (14.41) implies that they admit an explicit expansion 

*(v,0) = ^# n eW(v)e <w »*, (4.7) 

n 

where, we have suppressed the superscript + because from now on we will only work with 
positive frequency solutions. The physical inner product is given by 

(*i|tf 2 ) e = B(u)*iM ) * 2 Mo) (4-8) 

ve{±\e\+An;n&} 

for any <p . The action of the Dirac observables is independent of e, and has the same form 
as in the WDW theory: 

|%*(v,0) = e iVS |«| *(v,<f> ), and p^(u, 0) = -ih^^- . (4.9) 

The kinematical Hilbert space H^ 1 is non-separable but, because of super-selection, each 
physical sector 7ip hy is separable. Eigenvalues of the Dirac observable \v\^, constitute a 
discrete subset of the real line in each sector. The set of these eigenvalues in different 
sectors is distinct. Therefore which sector actually occurs is a question that can be in 
principle answered experimentally, provided one has access to microscopic measurements 
which can distinguish between values of the scale factor which differ by about a Planck 
length. This will not be feasible in the foreseeable future. Of greater practical interest are 
the coarse-grained measurements, where the coarse graining occurs at significantly greater 
scales. For these measurements, different sectors would be indistinguishable and one could 
work with any one. 



V. LOOP QUANTUM COSMOLOGY: NUMERICAL ISSUES 

As we saw in section [IV] physical states can be readily constructed from eigenfunctions 
of the difference operator 0. In the first part of this section, we study properties of these 
eigenfunctions. We show that admits only normalizable eigenfunctions with discrete eigen- 
values and numerically construct an orthonormal basis. In the second part we use this basis 
to construct and analyze physical semi-classical states. 



A. Spectrum of G 

Eigenfunctions of are solutions to the difference equation (14. ip . Since its coefficients 
are real, any eigenfunction can be expressed as a complex linear combination of real eigen- 
functions. Therefore, it will suffice to restrict ourselves to real eigenfunctions. 

Consider a generic lattice C e , i.e., a lattice where e does not equal or 2. Since C £ = 
C\ e \ U£-\ E \, to obtain an eigenfunction which is symmetric under orientation reversal v — > — v 
it suffices to solve the eigenvalue equation just on C\ e \ and then reflect it. Because is 
invariant under orientation reversal, the reflected function is automatically an eigenfunction 
with the same eigenvalue, but supported on C-\ e \. Furthermore, if the original eigenfunction 
is normalizable on C\ e \, the reflected one is normalizable on £_i e i whence the sum is a 
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symmetric, normalizable eigenfunction on C £ . Therefore, in what follows, for a generic e we 
will restrict ourselves to C\ E \ and examine all eigenf unctions. Since lattices with e = or 2 
are symmetric under reflection, on these exceptional lattices we will restrict ourselves only 
to symmetric eigenf unctions. 

Every such eigenfunction il) u (v) is uniquely determined by its 'initial' values ip u (e + 4n), 
t])u{£ + 4(n + 1)) for some integer n. 7 For later convenience, we note that the initial data 
can be represented by a pair of real parameters b G R and 9 6 [0, it] 

V> w (e + 4n) = bcos(6) , ip u (e + A(n + 1)) = bsm(6) . (JSHD 

Thus on generic lattices, for any u, the eigenspace is 2- dimensional and the degeneracy is 
parameterized by b, 9. On the exceptional lattices, the symmetry requirement imposes an 
additional constraint which determines ip u (e + 4) as function of ipu(e). Therefore on these 
lattices the eigenspace is only one dimensional. 



a) b) 




FIG. 2: Examples of normalizable and generic eigenfunctions of 0. The plot a) presents eigenfunc- 
tions supported on the lattice with e = 0.5. Since the scale on the y-axis is logarithmic, it is clear 
that in the asymptotic region generic eigenfunctions diverge exponentially while the normalizable 
ones decay exponentially. Plot b) shows eigenfunctions supported on the e = lattice. Since the 
scale on the y-axis is linear this plot brings out the detailed behavior of eigenfunctions in the region 
where \ipu\ is small. 

The generic behavior of eigenfunctions ipuj{v) is shown in Fig[2l One can distinguish three 
regions on the positive (or negative) v axis in each of which the eigenfunction ij) u (v) has a 
qualitatively different behavior. 

(i) Genuine quantum region: \v\ < Vb in which |^| grows/decays exponentially. For 
e = or 2, always grows with |i>|. Vb turns out to be proportional to uj. 



7 If e 7^ 2, the eigenvalue equation (|4.ip recursively determines ipu( v ) on the entire C+\ s \- If e = 2 the 
coefficients C~(2) and C + (— 2) in Eq. (|4.1[) vanish. Therefore, we can calculate ipu>{v) only on half of the 
v axis. However, values of 4'u>{' u ) on the other half are determined by the symmetry condition. 
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(ii) 'Semi-classical' region: t> fe < |u| < v r in which ip u oscillates. v r turns out to be 
proportional to uj 3 / 2 (and approximately equals the maximal \v\ of a classical universe 
of momentum p^ = ftw). 

(iii) Classically forbidden region: \v\ > v r , where ip u grows/decays exponentially as \v\ 
increases. 

(However, numerical simulations show that the distinction between first two regions gets 
blurred for eigenf unctions corresponding to uj < 5.) Note that for a generic lattice, the 
eigenfunction may decay in one asymptotic region (say, v — > oo) but grow in the other 
(v — > — oo). Eigenfunctions growing on either side fail to be normalizable, whereas the ones 
decaying for both signs of v are normalizable. While eigenfunctions exist for any uj, they 
decay on both sides only for certain discrete values of uj. Consequently the spectrum of 9 is 
discrete and the normalizable eigenfunctions form a set of zero measure in the space of all 
eigenfunctions. We will now describe the search algorithm to find them. We begin with the 
simpler case of e = or 2 lattices and then discuss the more subtle case of generic lattices. 

Fix e = or 2. Because of symmetry, it suffices to analyze the behavior of ip u (v) just 
for v > 0. Let us focus on an interval W = [a>i,a>2] of frequencies and fix a point v' G C £ 
such that v' ^> v r {uj 2 ). Eq. (14. II) implies that the value ip u (v') is a continuous function of 
i><j(v)\ v =£ and uj. One can fix the initial value ip u {e) to 1, thus leaving the dependence only 
on uj. Numerical inspection shows that ip^iv) changes sign quasi-periodically as uj increases 
(see the right plot of Fig. [2]). Let us take one of the (possibly many) values uj n>v i such that 
ifttj ,{v)\ v =v> — 0. The limits 

uj n = lim uj n ^> (5.2) 

are the only eigenvalues corresponding to normalizable eigenfunctions. In practice the values 
u n y for v' ~ 1.3v r approximate the limiting value u n with precision 1CT 16 . 
In actual calculations the following algorithm was applied: 

(i) First we consider a set of frequency values uji uniformly distributed within the interval 
[0, o; max ]. For each fixed v', the separation uj i+ i—Ui was chosen to be much smaller than 
separation between values of uj at which ip w {v') vanishes. (In practice the separations 
turned out to be greater than 1 for uj < 3 x 10 5 .) 

(ii) Whenever the change of sign between ip u>i+1 (v') and ^.(v') was detected, the value 
uj n y corresponding to the root of ^{v') was found via bisection method. 

For this scan, c^Wx was chosen as 3 x 10 5 . The first of the two steps in the algorithm ensures 
that all eigenfunctions in the chosen interval have been found. 

Let us now consider a generic lattice. Now there are two factors which complicate the 
task of finding normalizable eigenfunctions. First, the eigenspaces are two dimensional, 
parameterized by b,9 as in fl5.ll) . Therefore, even if we fix the normalization freedom by 
setting 6 = 1, for each frequency uj we have a 1 -parameter family of eigenfunctions, labeled 
by 9 (rather than a single eigenfunction as on e = or 2 cases). Secondly, since the desired 
eigenfunctions do not have to be symmetric on C\ £ \, now we have to analyze the behavior 
of ipui when |u| ^> v r separately for positive and negative values of v. We therefore modified 
the algorithm specified for the e = or 2 lattices as follows. 

First, keeping uj fixed and varying 9 instead of uj in the above procedure, we searched 
for ipuj{ v ) which decays on the negative v side. For this we probed the domain of 9 in 100 
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FIG. 3: a) Distribution of lowest (a; < 50) eigenvalues (corresponding to normalizable eigen- 
functions) for different e. Eigenvalues are divided into three classes according to the behavior of 
eigenf unctions, which are i) suppressed for v < 0, ii) suppressed for v > 0, and hi) large only near 
v = 0. b) Rescaled spectral density 2it/(oj n — u> n _i) (i) is compared with periods T${u} n ) in cf> (iii) 
of expectation value of |«L of the coherent states evolved from the WDW initial data (see section 
IV Bp . The solid line (ii) shows the large u limit of the rescaled spectral density. 



points uniformly distributed within [0, ir] and, whenever sign of il> u (v') changed, we narrowed 
the choices of 9 using bisection. The analysis shows that for each uj there exists unique 9 
such that ij) w satisfies this condition. With 9 so determined, we again have a 1-dimensional 
eigenspace ^(v) for each uj. Therefore, we could now apply the procedure used for e = or 
2 to look for eigenfunctions which decay on the positive side. Thus the two complications so 
to say compensate one another: while there is an additional, 1-dimensional freedom (labeled 
by 9) in the choice of eigenfunctions, since the requirement that the eigenfunctions decay on 
the negative v side is now decoupled from the requirement that they decay on the positive 
side, we have an additional constraint (which determines 9 for any given uj). Therefore, the 
only modification to the algorithm used in the e = or 2 cases was to first determine 9. 

To gain insight into the qualitative features of eigenfunctions and eigenvalues, this method 
was first applied to uj G [0, 50] by choosing C + \ £ \. Normalized eigenfunctions decay exponen- 
tially as \v\ increases for both positive and negative v. It was found that these eigenfunctions 
generically have an additional feature: in the 'semi-classical' region described above, their 
values on the positive or negative sides of the v axis are suppressed relative to their values 
on the other half, the difference in the amplitudes growing exponentially with uj. Since we 
will be primarily interested in large d, if the initial data are specified on the side where ip w 
is large, the suppression on the other side enhances the numerical errors making the results 
unreliable. For these cases, the search of normalizable eigenfunctions was performed again, 
now starting from the side where ip^ is small. 

After implementing all these precautions to control errors, a much more exhaustive search 
for eigenvalues and eigenfunctions was made. The results can be summarized as follows. 

(i) The spectrum of is discrete for all lattices. Each of its eigenvalues is non-degenerate. 
Normalizable eigenfunctions decay exponentially as |i>| tends to infinity. 
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(ii) For generic e, the restrictions to sublattices £±|e| of eigenf unctions if) n (v) are in general 
strongly suppressed for one sign of v. The side on which suppression occurs depends 
on parity of n. 

(iii) The density of eigenvalues has following features. As shown in Fig. [3^ it is roughly 
independent of e except that for generic lattices it is twice as large as that for excep- 
tional ones. It slowly grows with uj (see Fig. [3b). The best fit to the data in Fig. [3b 
leads to the behavior in large u limit as: 

u n — ~~ ► — -, — i 7-57 (5-3) 

a In \u n /p\ 

with a = [0.0259272 ± (5 x KT 7 )]^ 1 / 2 and (3 = [0.04412 ± (1 x IQ-^G 1 ' 2 . 

(iv) A state ty(v) which is sharply peaked at some uj* can be well- approximated by a 
linear combination only of eigenfunctions with eigenvalues lying in a small compact 
interval around u*. If u* is very large the estimate given in (iii) above implies that 
the distribution of uj n is approximately uniform. Consequently, the wave function will 
be approximately periodic in <f> with a period 

T^K) « . (5.4) 

For an LQC physical state obtained from evolution of the initial data at <fi = <p 
corresponding to a WDW coherent state (see section IVBl) . the period of expectation 
value (\v\cj,) turns out to be in good agreement with T^. This provides an independent 
check on our numerical analysis. 



B. Evaluation of Semi-classical State and Observables 

The numerical method presented in previous subsection allow us to find all the eigenstates 
of O which span the physical Hilbert space. To obtain the normalized eigenbasis (which 
is symmetric under orientation reversal) we first note that in the expression of the norm it 
suffices to evaluate just a finite sum 

HnWl = Where N> J- ( 5 - 5 ) 

v e{±e+n; ~N<n<N} 

3 

since ip n {v) decay exponentially for \v\ > v r oc oji . With the basis at hand, one can construct 
physical states which are semiclassical at late times. As in k = case 0], this can be done 
in two different ways: 

(i) Direct evaluation of the wave function using (14.71) . and 

(ii) Evolution in of the initial data specified at = O using Eq. (I2.3ip . (The initial data 
can be chosen to be the same as that of a semi-classical solution to the WDW equation 
at some late 'time' (see section ITTTT) ). 

The direct evaluation of the integral expression (14. 7p already provides the full LQC solu- 
tion. However, we also used this expression to obtain just the initial data at (f> = (fi and 
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then evolved this data using (I2.3ip . Agreement between the two solutions provides an in- 
dependent check on numerics. Finally, we also evolved the initial data extracted from a 
WDW coherent state and used the resulting LQC solution thereby obtaining an indepen- 
dent check on completeness of the eigenbasis constructed in section IV Al The rest of this 
section provides the relevant details of the numerical implementation of this procedure. 

In order to evaluate the integral solution, we need the spectral profile \P n used in the 
expansion (I4.4p . Since we are primarily interested in states which are semiclassical at late 
times representing a macroscopic universe, we chose a Gaussian profile peaked around large 
uj*\ \P n = exp(— (uj n — u*) 2 / (2a 2 )) . Since the contribution from eigenf unctions corresponding 
to uj sufficiently far away from u* can be neglected, to calculate ^(v,<f)) via (14. 7p one has 
to only sum a finite number of terms. In numerical simulations the summation in (14. 7[) was 
restricted to n such that u;* — 10a < uj n < u* + 10a. 

To calculate the initial data from the WDW coherent state we evaluated (13.151) . Now 
the eigenf unctions tk(v) are appropriately normalized Bessel functions /C^. They were cal- 



culated using methods of Gil, Segura and Temme [24J. As in the initial data construction 
discussed in the last para, it is sufficient to restrict the domain of integration to a compact 
set [k* — 10a; k* + 10a]. Resulting integral was then numerically evaluated using a simple 
trapezoid method with set of 10 4 'probing' points distributed uniformly. On the one hand, 
this calculation immediately leads to the solution to the WDW equation presented in Fig. 
[TJ On the other, it provides us LQC initial data for any lattice C £ which was evolved using 

(USD- 

Eq. (12.311) constitutes a countable number of ordinary differential equations. Its domain 
in v is the lattice C e . However, as noted in section IV Al the symmetry of the wave function 
allows us to restrict the calculations to the sublattice £+\ £ \ for generic e and to the part 
v > in the exceptional cases e = or 2. Due to technical limitations the size of the domain 
was restricted by requiring that its elements i>, satisfy the inequality \v{ — e\ < 4iV ^> v r (u*). 
To ensure that the system remains closed, we impose boundary conditions on the outermost 
points — e\ = 4N in the generic case and at the right outermost point v , — e = 4N in the 
exceptional cases. Since these points lie deep in the classically forbidden region for all the 
eigenvalues contributing significantly to the state, one can safely impose reflective boundary 
conditions 

V{±4N + e,<p) = fytf(±4iV + e,0) = 0. (5.6) 

In actual simulations we chose AN pa 1.3v r (u*) which was sufficient for the boundary condi- 
tion to not affect the dynamics. 

The resulting finite set of equations was integrated using adaptive 4th order Runge- 
Kutta method. To estimate the numerical errors due to discretization in 0, restrictions 
calculated for different step sizes were compared using sup-norm 

||/||(0) = sup f(v u <P) . (5.7) 

\vi-e\<4N 

The step sizes A0 were chosen to satisfy the inequality 

||*A*-*A«/2|| < e ||*A*/2|| A0 (5.8) 

for small preset e, where \I/a</> and ^a</>/2 are profiles calculated with step sizes A</> and A0/2 
respectively. The dependence of H^a^ — ^A^H on the number of steps of the integration 
for one of the calculations is presented in Fig. HI It shows that the numerical errors mani- 



24 



fest themselves mainly in phases. The differences between the absolute values of the wave 
function profiles are approximately one order of magnitude smaller. Thus the expectation 
values and dispersions of observables |{*|^ are determined with much better precision than 
* itself. 




10 4 10 5 10 6 10 7 

M 



FIG. 4: Error functions H^m) ~~ ^ ll/ll^ll (upper curve) and IH^m)! ~~ \^\ ll/ll^ll (lower curve) 
are plotted as a function of the number of time steps. Here ^(m) refers to final profile of wave 
function for a simulation with M time steps, VP is the limit of final profile as 1/M — > calculated 
via polynomial extrapolation. In both cases, the evolution began at (ft = and the final profile is 
evaluated at (ft = —1. Values of parameters are p^ = 1000, Ap^/p^ = 0.018, v* = 4000 and e = 0. 

The resulting wave functions ^f(v, (ft) were finally used to calculate the expectation values 
{pct>), (H(/>) of observables defined by (14. 9p . With the inner product (\l/|\l/) e given by (14. 8 p 
they are equal to the following sums over C e ^ := {v = ±£ + 4n; —N < n < N} 

(#| \v\^\y} = (^i^); 1 B ( v )\ v \\^( v ^)\ 2 ( 5 - 9a ) 

(*|^|*> = (^l^); 1 B(v)*(v,<l>)(-iK)d4p{v,<j>) (5.9b) 

The dispersions corresponding to considered observables were calculated via 

(A^> 2 = (pi) - (P4,) 2 (5.10a) 
(A|%> = @J> - (H,) 2 (5.10b) 

where the expectation values (p|) and are of the form (15.91) . 

In actual simulations to* ranged from 500 to 10 4 , with relative uncertainties Aw/w between 
7 x 10~ 3 and 2.5 x 10~ 2 . Wave functions were calculated on lattices C e corresponding to 5 
different values of e. Peaks in v of initial profiles were chosen to be no smaller than half of 
the maximal v predicted by classical theory. For example for = 500, v* ~ 2000 whereas 
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FIG. 5: The absolute value of the wavefunction obtained by numerical evolution. For visual 
clarity, only values 1 1 QE r | > are shown. Parameters of the initial data are: pj = 5 x 10 3 , 

A-P4>/p<f> = 0.018, v* = 5 x 10 4 and e = 0. Because of bounces the quantum universe exhibits a 
'cyclic' character. 



for = 5000 the value of v* ranged between 5 — 6 x 10 4 . Evolution was performed in both 
forward and backward direction. Representative results of numerical evaluation of 
and the expectation values of \v\^ are shown in Figs. O and El Detailed discussion of the 
properties of \^(v , <f>)\, the comparison of classical and quantum evolution and a summary 
of our results is presented in Sec. IVIBI 



VI. PHYSICAL IMPLICATIONS 



In this section we discuss the physics of our numerical simulations. We first obtain the 
quantum corrected effective Friedmann equation and show that it is an excellent approxima- 
tion to the numerical quantum evolution. We then list the numerical results and compare 
and contrast the exact quantum evolution, the effective theory and the classical Friedmann 
dynamics. 



A. Effective Equations 

The right hand side of the standard Friedmann equation 



B» S (5)> S (6.1) 
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FIG. 6: Expectation values and dispersion of for the wavefunction in Fig. [5] are compared 
with the classical trajectory and the trajectory from effective Friedmann dynamics obtained from 
Eqs. (l6.3p and (|6.4j) . The classical trajectory deviates significantly from the quantum evolution at 
Planck scale and evolves into singularities. The effective trajectory provides an excellent approxi- 
mation to quantum evolution at all scales. 



vanishes if and only if p = 3/8irGa 2 . In any classical solution, at this point the scale factor 
reaches its maximum value, the matter density its minimum value, and the classical dynamics 
exhibits a turning point from an expanding phase to a contracting phase. Our numerical 
simulations show that quantum dynamics retains this turning point but also gives rise to 
additional turning points, resolving the big bang and the big crunch singularities (see Fig. 
|5]) . To gain an analytical understanding of the physics underlying this phenomenon, in this 
subsection we will obtain the leading LQC corrections to the classical Friedmann equation. 
This quantum corrected or 'effective' Friedmann equation already suffices to describe the 
behavior of the peak of wave functions that we found numerically in section |V] (see Fig. [6]). 

The procedure for obtaining these effective equations is the same as that in the k=0 
case. Let us begin with a brief summary of the main ideas. One begins with a geometrical 
formulation of quantum mechanics in which the space of pure quantum states is represented 
an infinite dimensional symplectic manifold r quan — called the quantum phase space (see, 



e.g., (25]). The quantum phase space has the structure of a fiber bundle: the base space 



is the classical phase space r c i ass and the infinite dimensional fiber over any point (q Q , p Q ) 
of r c i ass is the space of quantum states in which the expectation values of the canonically 
conjugate operators (q, p) are (q ,Po)- 8 Interestingly, the exact quantum dynamics provides 
a Hamiltonian flow on the symplectic manifold r quan (the corresponding Hamiltonian being 



In LQC (and LQG) an important subtlety arises because there is no operator corresponding to the 
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just the expectation value function of the quantum Hamiltonian operator). To obtain the 



desired, first order quantum corrections, one finds a cross-section of T 
ding ofT c i ass into T, 



quan 



-i.e., an embed- 



sense 



quan to which this flow is approximately tangential in a well-defined 
This approximate quantum Hamiltonian flow unambiguously projects down to 
r c i ass and provides the desired corrections to classical equations of motion. In a certain sense 
this procedure encapsulates the more familiar 'effective action' calculations in the Hamil- 
tonian framework J271I . In LQC, these quantum corrections have been obtained for various 
matter sources [21 



281 ] . where judiciously chosen generalized coherent states are used to de- 



fine the required embedding of r class into r quan . Just as the standard effective action refers 
to the in and out vacuum states, the final effective Hamiltonian (and thus the Friedmann 
and Raychaudhuri equations) depend on the specific choice of coherent states. However, the 
first order corrections we are interested in are insensitive to these details. In our case, the 
resulting effective Hamiltonian is: 



H 



a 



eil' 



cff 



A(v) 



16nG 16nG 



sin ^(c--)-sm (—) + (! + r) — 



+ 



6 J K J 2 
(6.2) 

where A(v ) denotes eigenvalues (I2.26|) of A . 

Since the Friedmann equations involve a/ a, to obtain modifications, we first derive the 
Hamilton's equations of motion: 



and 



v = {v, H eS } 



8n"fG dHes dv 
3 dc dp 



.T^M f^)- 1 K ^ \ v \ V3 sin ^ _ |) cos A(c _ |) 



{(j), H eS } = { — - — J B(v)p^. 



(6.3) 
(6.4) 



Using the constraint equation 7i cff = 0, i.e. 



sin 



,2 - , 



fi(c-%) = Uin' (i^)-(l +7 2 ) 



(6.5) 



we can eliminate the dependence on the connection c in (16.31) and obtain the desired 
quantum-corrected Friedmann equation 



H 2 



Yfi 2 A(v) 

(87r 7 £; 
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8ttG 



'2B(v) 



A(v) 



(6.6) 



which, as with the classical Friedmann equation (16. ip . involves only v, v and p^. To take the 



configuration variables, i.e., the connections. One has to use holonomies instead. This issue is handled in 
26|. 
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classical limit, we note that the volume is given by (87T7/6) 3 / 2 (\v\/K) £|j and that (the area 
gap and hence) p goes to zero in this limit. Therefore, in the limit (16.61) reduces precisely 
to (16.11) . Terms containing p represent the quantum geometry corrections. 

In view of the fact that the quantum equations are invariant under the orientation reversal 
map H^{y) = v), it will suffice to restrict ourselves to v > 0. Now, for v > 1 we have 



m = -^m (6.7) 



and for v 3> 1, we have: 



1/3 

73/2^2 \K 



B(v) = — + 0(v- 3 ) (6.8) 

v 



Therefore, the leading order quantum corrected equation is given by: 




X 



l-siu 2 (^)-(l + 7 2 )^-^) + 0(,"») (CO) 

z Pcrit , 



where p denotes the eigenvalues of the energy density operator p = p|/|p| 3 and as in the k=0 
case we have set p crit = 3/(16TT 2r y 3 G 2 h) ~ 0.82p P i. Although for brevity we have kept the 
sine functions in this equation, since v _1 ~ p 3 , to the leading order considered here one only 
needs to keep terms 0(p s ) in their Taylor expansions. Finally, note that on substituting 
£ = in (I6.9p . we immediately obtain the effective Friedmann equation of j3| for the k=0 
model. 

To probe the possible turning points, it is useful to rewrite the modified Friedmann 
equation (I6.9P in the following form 

H 2 = ^ (p - Pi) (—{pu - p)) + 0(v- 3 ) (6.10) 



Pcrit 



with, 



8ttG^ 2 p 2 V 6 

and 



P2M = Palt (l - Sin 2 (^) - (1 + 7 2 )^) • (6.12) 



Note that while p cr it is a constant, pi,p2 are functions of v (and p is a function of v, 
Along each dynamical solution of the effective evolution equations (16.31) and (16.41) . p, p 1; p 2 
all evolve. In the effective dynamics, along any given dynamical trajectory turning points 
occur when p = p\ or p = p 2 . Plots of solutions to the effective equations show that the 
classical recollapse occurs when p — pi, where the universe reaches its maximum radius a max 
and minimum density p m i n , and the quantum bounce occurs when p = p 2 where the universe 
reaches its minimum radius a m j n and maximum density p max - Now, an examination of the 
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expressions of p\ and p 2 show 
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Pcrit 




Numerical simulations showed that our notion of 'semi-classicality at late times' is surpris- 
ingly weak. For example, in the simulation with p^ pa 5 x 10 3 h~ (in the classical units c=G=l), 
the universe grows only to a maximum radius of pa 23£p\ before undergoing the classical rec- 
ollapse. Even for this small universe, effective equations predict that the density p m \ n at 
the recollapse should agree with the classical Friedmann formula p m i n = 3/87rGa max to one 
part in 10~ 5 and the density p max at the quantum bounce would equal the critical density 
Pcrit — the density at the bounce in the k=0 models — to within a couple of percent. These 
predictions are borne out in the numerical simulations of the exact LQC equations. Since 
a max scales as (p^) 1 ^ 2 and a min as (p</,) 1//3 , effective equations imply that both approximations 
improve as p^ increases and become almost exact for universes which grow to interesting 
macroscopic sizes. 

In particular, as one would hope, the effective theory accurately reproduces the predic- 
tions of classical general relativity in the low curvature regime. Yet, already the leading 
order correction from quantum geometry is strong enough to resolve singularities and re- 
place them with a bounce. It is also noteworthy that for universes which grow to large 
macroscopic sizes, the density at the quantum bounce is universal, and equals that in the 
k=0 models. These general features represent key predictions of the numerical evolution of 
the full quantum equations. The effective theory provides a physical understanding of how 
they come about. In particular it shows that the most important corrections come from 
the gravitational part of the quantum Hamiltonian constraint. Indeed, since we need only 
the leading, classical part of B(v) in this analysis, the quantum modifications of the matter 
Hamiltonian play no role in the leading order corrections discussed here. 

In the last part of our above discussion we have ignored terms 0(v~ 3 ). A priori it is pos- 
sible that the sum of these terms is not negligible and could even dominate the leading term. 
However, a comparison with numerical simulations shows that the predictions obtained from 
just the leading order quantum corrections in the effective theory accurately describe the 
full quantum dynamics of states which are semi-classical at late times and the accuracy im- 
proves as p<p increases. Thus, although one could not draw definitive conclusions just from 
effective equations, when used in conjunction with numerical simulations, they provide an 
easily manageable and powerful tool to probe quantum geometry effects for universes which 
grow to a max ~ 25£pi or more. 



Main results on quantum dynamics can be summarized as follows. 

1. Consider a classical solution which evolves from the big-bang to the big crunch, reach- 
ing a large maximum radius a max . Fix a point on this trajectory where the universe 
has reached macroscopic size and consider a semi-classical state peaked at this point 
(see section |VB|) . Such states remain sharply peaked throughout the given 'cycle', i.e., 
from the quantum bounce near the classical big-bang to the quantum bounce near the 
classical big-crunch. The notion of semi-classicality used here is rather weak: these 



B. Results 
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FIG. 7: a) Zoom on the portion near the bounce point of comparison between the expectation values 
and dispersion of v\s, the classical trajectory and the solution to effective dynamics. At large values 
of \v\fj) the classical trajectory approaches the quantum evolution. Values of parameters are the 
same as in Fig. [6l b) The behavior of expectation values of p for different values of p*^ are shown. 
On each solution, these are bounded between a p max and a p m m- For a universe peaked at large 
values of p^, p max « p crit . 



results hold even for universes with a max pa 25£p\ and the 'sharply peaked' property 
improves as a max grows. 

2. The trajectory defined by the expectation values of the Dirac observable v\^, in the 
full quantum theory is in good agreement with the trajectory defined by the classical 
Friedmann dynamics until the energy density of the scalar field becomes comparable 
to the maximum energy density p max ~ p cri t ~ 0.82p P1 . Then the classical trajectory 
deviates from the quantum evolution. In the classical solution, the matter energy 
density keeps increasing on further evolution, eventually leading to a big bang (re- 
spectively, big crunch) singularity in the backward (respectively, forward) evolution, 
when v — > 0. The situation is very different with quantum evolution. Now the uni- 
verse bounces at p = p max , avoiding the past (or the big bang) and future (or the big 
crunch) singularities. 

3. The expectation values and relative dispersions of remain constant during different 
stages of evolution. Thus the expanding and contracting branches correspond to the 
same value of (p^). Further, as a check on numerics we verified that the norm of the 
states is also preserved during the entire evolution. 

4. After the quantum bounce the energy density of the universe decreases and, when 
P "C Pmax, the quantum evolution is well-approximated by the classical trajectory. On 
subsequent evolution, the universe recollapses both in classical and quantum theory 
at the value v = f max when energy density reaches a minimum value p m i n . 

5. The trajectory obtained from effective Friedmann dynamics (I6.9P is in excellent agree- 
ment with quantum dynamics throughout the evolution. (See Figs. [6] and [7^.) In 
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particular, the maximum and minimum energy densities predicted by the effective 
description agree with the corresponding expectation values of the density operator 

p = p^/\p\ 3 computed numerically Evolution of the expectation values of p with \v\,p 
is shown in Fig. [7b. 

6. For quantum states under discussion, the density p max is well approximated by p crit ~ 
0.82ppi up to terms 0(£^/a ^^ ), independently of the details of the state and values 
of Pip. (For a universe with maximum radius of a megaparsec, ^pi/a^in ~ 10 -76 .) The 
density p min at the recollapse point also agrees with the value (3/(87rGa^ ax ) predicted 
by the classical evolution to terms of the order 0(£ P1 /a^ ax ). Furthermore the scale 
factor a max at which recollapse occurs in the quantum theory agrees to a very good 
precision with the one predicted by the classical dynamics. 

7. The relative dispersion of increases — but very slowly — as one evolves through 
many cycles. Effective Friedmann dynamics provides insight into this behavior of the 
quantum state. Let us first consider two nearby solutions to the effective equation 
with slightly different p^ but with same value of v at a chosen = O . Then the 
relative difference between v of the two solutions after one cycle can be estimated 
using Eqs. (16.31) and (16. 4p as 

— « V^G a ^ (6.14) 

v P0 

where Sp^ is the difference between values of p^ of the two effective trajectories and 
y/G a = 0.0259272 ± (5 X 10~ 7 ) (see 15.31) . This estimate was found to provide a good 
upper bound on the relative dispersions computed using numerical evolution of the 
quantum state. 

8. The state remains sharply peaked for a very large number of 'cycles 7 . This number 
can be estimated using Eq. (l6.14l) . Consider the example of a semi-classical state with 
an almost equal relative dispersion in p^ and \v\cf, and peaked at a large classical 
universe of the size of a megaparsec. When evolved, it remains sharply peaked with 
relative dispersion in \v\$ of the order of 10~ 6 even after 10 50 cycles of contraction 
and expansion! Any given quantum state eventually ceases to be sharply peaked in 
\v\tj, (although it continues to be sharply peaked in p^). Nonetheless, the quantum 
evolution continues to be deterministic and well-defined for infinite cycles, i.e., on the 
entire real line of the emergent time 0. This is in sharp contrast with the classical 
theory where the equations break down at singularities and there is no deterministic 
evolution from one cycle to the next. In this sense, in LQC the k=l universe is 
cyclic, devoid of singularities. As in the k=0 case, this non-singular evolution holds 
for all states, not just the ones which are semi-classical at late times. There is no fine 
tuning of initial conditions. Also, there is no violation of energy conditions. Indeed, as 
discussed in section IVI A[ quantum corrections to the matter Hamiltonian do not play 
any role in the resolution of the singularity. The standard singularity theorems are 
evaded because the geometrical side of the classical Einstein's equations is modified 
by the quantum geometry corrections of LQC. 



9. In the k=l model, certain effective equations have been written down and used to 
predict non-singular bounces 29] in the broad framework of LQC. The presence of 
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these bounces was also used to study the onset of a successful period of inflation 
in closed models j30|. How do these analysis compare with that presented in this 
paper? There are two important differences. First, these works focused on the matter 
part of the Hamiltonian constraint and made a crucial use of quantum corrections to 
the matter Hamiltonian arising from the use of representations (of SU(2)-holonomies) 
labeled by large values of j. Second, these large j representations were used only 
in the matter part of the Hamiltonian and not in the gravitational part. Since then 
Perez [l7[ has shown that mathematical consistency requires us to use the j = 1/2 
representation in 3-dimensional gravity. He also argues that the same should hold 
in 4 dimensions. In any case, while the use of, say j = 1 representation could be 
justified because the model has no spinor fields, the use of large j values appears to 
be unnatural and needs an independent justification which is still lacking. Similarly, 
without an independent justification, an asymmetric treatment of the gravitational 
and matter appears to be ad- hoc [16( , somewhat similar to using two different metrics 
on the right and left hand side of classical Einstein's equations. Our analysis used the 
fundamental, j = 1/2 representation for geometry as well as matter. Our numerical 
simulations as well as effective equations show that with this choice the modifications 
of the matter Hamiltonian play no role. This is consistent with findings in the older 
literature; indeed this is the reason why large j values were used there. In our analysis, 
it is the modification of the gravitational part of the Hamiltonian constraint that plays 
a key role in the singularity resolution. In the older literature, by contrast, these 
modifications were ignored. Nonetheless, basic physical ideas in these older works are 
intriguing and it would be interesting to reanalyze those issues using the Hamiltonian 
constraint introduced in this paper. 

We will conclude this section by clarifying two issues that have arisen from a recent work 
[11] addressed to computational physicists, particularly numerical relativists. Although this 
discussion refers only to the k=0 model considered in [H, 0, [3|, it is included here because 
the same issues can arise in the k=l model. 

The phrase "These bounces can be understood as spurious reflections" in the abstract 



(and again in the body) of [llj was interpreted by some as suggesting that the bounces 
reported in [0, @, 3 were artifacts of bad numerics. This is certainly not the case: Not 
only were those simulations performed with all the due care but our result that the LQC 
equations predict a genuine, physical bounce was in fact reproduced in the first half of 
[ll| . Indeed, this part is a nice summary of our numerical results geared to computational 
physicists. From discussions with the author we understand that the intent of that phrase 
was to say: 'had the physical problem been to solve a wave equation in the continuum 
and had one used non-uniform grids, one would also have found bounces which, from the 
perspective of continuum physics of this hypothetical problem, would be interpreted as 
spurious reflections in finite difference discretizations'. This is likely to be an illuminating 
point for computational physicists but is not physically relevant in LQC where the basic 
equation is a difference equation. 

The second part of the paper considers a modification of the quantum Hamiltonian con- 
straint by "adding ad-hoc higher order terms." It is then suggested that such modifications 
could remove the bounce. Let us analyze the issue from a mathematical perspective even 
though the analysis and conclusions have no obvious physical significance since the modifi- 
cations do not result from any systematic, physical considerations. Then, since the physical 
state is symmetric under orientation reversal, simulations reported in ll[ imply that the 
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bounce would not disappear but change its character. In the physical solution there would 
again be a pre-big-bang, contracting branch which would be joined in a deterministic fashion 
to a post-big-bang expanding branch. However, now the two branches will meet at v = 0. 
Although the ensuing differences are not trivial, the qualitative picture is not changed even 
by these ad-hoc modifications. 9 

VII. DISCUSSION 
A. Key Features 

Key features of the k=l model can be summarized as follows. 

i) The scalar field serves as emergent time at all three levels: classical general relativity, 
WDW theory and LQC. In the classical theory, every solution undergoes a recollapse but 
<j) remains single valued. Each solution begins with a big bang and ends with a big crunch 
and <p ranges over the entire real line irrespective of the constant of motion p^. In the 
WDW theory and LQC the form of the Hamiltonian constraint operator implies that can 
serve as emergent time also in the quantum theory. Situation with the range of <j) i n the 
WDW theory is the same as that in classical general relativity since the singularities are 
not resolved, i.e., since quantum dynamics can not unambiguously evolve the state across 
these singularities. In LQC on the other hand the singularities are resolved and the quantum 
evolution across the putative classical singularities is deterministic. The range of <ft continues 
to be the entire real line. 

ii) In LQC, three sets of results show that the big bang and the big-crunch singularities 
are both resolved. First, the LQC effective equations do not break down. Rather, while the 
classical Friedmann equation has only one root at which H 2 = (a/ a) 2 vanishes, quantum 
corrections introduce a second root when the matter density enters the Planck regime, 
altering classical dynamics and giving rise to bounces. The second result refers to full 
quantum dynamics: in contrast to the WDW theory, in LQC every state in 7i p h y has a 
well-defined unitary evolution for the full range (— oo, oo) of the 'emergent time' <fi. The 
third set of results is more detailed. It again involves the full dynamics of LQC but only 
semi-classical states in 7i p h y . Consider a classical trajectory in the v-<p plane in which the 
universe evolves to a macroscopic maximum size. The classical Friedmann equation implies 
that the universe attains its maximum volume V mgx at the recollapse point and this value is 
related to the constant of motion p^ via 10 

VW = (16rrG/30M « 0.6pJ. (7.1) 

Hence, in the solution under consideration p^ has to be large. Consider a point on this 
classical trajectory at a late time when the volume is macroscopic and a semi-classical state 



9 Indeed, as the numerical simulations of [2|,|3( show, the 'reflection' and transmission' phenomenon discussed 
in the second part of occurs even in standard LQC. It does not have deep physical significance in the 
final picture because physical states are symmetric under the orientation reversal map II. 

Here and in what follows, numerical values are given in the classical units G=c=l. In these units has 
the same physical dimensions as h and the numerical value of H is 2.5 x 10 _66 cm 2 . 
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in which Dirac observables and \v\^, are peaked at this point at 'time' = O . When 
evolved, this quantum state remains semi-classical, sharply peaked at the classical solution 
under consideration during the entire cycle, except near the big-bang and the big-crunch. 
There, because of repulsive effects of quantum geometry, the wave function bounces when 
the peak reaches a minimum volume: 

Knin = { AlXG ^ f 2 P^ « (1-28 x lfr 33 cm) P4> (7.2) 

Thus, the wave function does not enter a neighborhood of the classical singularity. The size 
of this neighborhood is dictated by the value of the constant of motion p^ and can be much 
larger than the Planck size. Finally, note that in the mathematical limit in which the area 
gap A is taken to zero, V m { n vanishes. More generally, this is also the limit in which LQC 
reduces to the WDW theory. Thus, one can use A as a knob to turn on or off quantum 
geometry effects and understand the 'mechanism' behind the singularity resolution. 

iii) To better understand the physical conditions at the quantum bounce, let us consider 
a couple of examples. Consider first a quantum state describing a universe which attains 
a maximum radius of a megaparsec. Then the quantum bounce occurs when the volume 
reaches the value Vmi n ~ 5.7 x 10 16 cm 3 , some 10 115 times the Planck volume. As a second 
illustration, consider a semi-classical state representing a large universe whose density at 
the classical re-collapse is about the current density of our universe p ~ 9.7 x lCT 30 gm/cc. 
Then, the quantum bounce would occur at V m i n ~ 1.4 x 10 24 cm 3 . Thus the 'quantum' 
or the 'Planck' regime is defined not by the volume of the universe but by the value of 
the matter density, or space-time curvature. In universes which grow to large macroscopic 
sizes, these quantities can attain Planck scale even when the volume of the universe is large. 
Figures quoted above were arrived at using our model where the scalar field is massless. The 
presence of potentials could significantly modify their values. Still, since the values of Vmm 
are so huge, these considerations are useful in drawing qualitative conclusions. For example, 
they suggest that in universes which grow to macroscopic sizes, the so called 'dj-effects' 
associated with modifications of the matter Hamiltonians due to quantum geometry will 
not be dynamically significant in homogeneous models (unless one considers astronomically 
large — and hence implausible — values of j). 

iv) Since detailed predictions were obtained only for states which are semi-classical at 
late times, it is interesting to ask how semi-classical these states have to be. How quantum 
mechanical can we make the parameters of the universe, still keeping the quantum state semi- 
classical in the sense used in this paper? The typical values of p^ used in the simulations 
reported in section |V] was 5 x lO 3 ^. These universes evolve only to a maximum volume of 
Vmax ~ 2.3 x 10 5 £p 1 before undergoing a recollapse. Results of our numerical simulations show 
that the necessary semi-classical considerations hold already for such small universes. In 
particular, the maximum value p max of the matter density is well approximated by the critical 
value p cr it = 0.82 pp\ already for these universes. A combination of numerics and effective- 
equations shows that the approximation becomes increasingly better as one considers larger 
and larger universes. Similarly, our numerical results have shown that the recollapse occurs 
at the classically predicted values already for these universes. In these cases, the matter 
density even at the recollapse point is quite high, approximately 2.2 x 10~ 4 ppi. Thus, 
several interesting phenomena occur in a rather small interval (of just four to five orders 
of magnitude) of density and volume. As our detailed plots show, there is a very narrow 
range of these parameters in which quantum geometry effects become significant. They 
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grow extremely quickly, overwhelm the classical attractive force, cause the bounce and then 
become insignificant very quickly again. 

v) We would like to emphasize that while quantum geometry effects resolve classical sin- 
gularities, we did not predict the emergence of classicality at late times within any given 
cycle. 11 These two issues are logically distinct. Indeed, in our detailed analysis we simply 
restricted ourselves to a single cycle and to states which are semi-classical at late times 
therein. At first, it may appear that the problem of actually specifying such states would 
be impossibly difficult in more realistic models since such specification would have to incor- 
porate all the complexities that have developed during the epoch during which the universe 
grew to a macroscopic size. However, as discussed above, this 'macroscopic size' can be very 
small and the required specification of a semi-classical state could be done at a relatively 
early time before complicated structures develop. For example, in an inflationary scenario 
one could specify the state immediately after the end of inflation, or perhaps even before the 
onset of inflation since, in the current observationally favored scenarios, matter density is 
significantly smaller than the Planck density even at the onset. Thus, while the conceptual 
issue of singling out a preferred family of states using general principles remains largely 
unexplored, there do not appear to be any 'practical' difficulties in specifying the state. 

vi) Since both the big bang and the big crunch singularities are resolved, the quantum 
wave function evolves through infinitely many classical cycles. Thus, in the k—1 model, the 
quantum space-time of LQG is vastly larger than that the classical space-time of general 
relativity. The issue of emergence of semi-classicality raised above can now be elevated to 
the infinite history of the quantum universe. Let us then begin a classical solution in which 
the maximum volume of the universe is large, say more than a megaparsec, and consider a 
quantum state which is peaked at a point on this classical trajectory at a late 'time', <fi = <p . 
As discussed in section [VT] this state has interesting properties. Except near the classical 
singularities, it will remain sharply peaked not only on the given classical solution during 
its cycle but also for over 10 50 cycles resulting from quantum bounces. However, because 
the eigenvalues of — id^ = VQ are not exactly evenly spaced, eventually the wave function 
will spread and cease to be semi-classical. The issue of converse is intriguing. Consider 
any state which is sharply peaked at a large value of the constant of motion p^ but has a 
large spread for the volume operator \v\<p at the initial instant of 'time' <fi = <p . Initially, 
such a state is not semi-classical. However, would the LQC dynamics evolve it to a state 
which is eventually sharply peaked at a classical trajectory? The answer appears to be in 
the affirmative [31] . If so, all states which are peaked at a very large value of p^ would 
eventually become semi-classical. In this precise sense, semi-classicality would be generic. 



B. Classical recollapse from LQC 

As discussed in section [J a major challenge to any background independent quantum 
gravity approach, such as LQG, is to ensure that there is a sufficiently large semi-classical 
sector. LQC offers a non-trivial context to probe this issue. Are the quantum geometry 
effects subtle enough to dominate near classical singularities but turn themselves off on 



As discussed in the last section of [2|, the bounce picture has a suggestive, intriguing relation to the 
Hartle-Hawking [23| proposal for the wave function of the universe. It may well suggest an avenue to 
address this issue. 
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large scales? As results of [2] show, this is a delicate issue. Indeed, the evolution generated 
by the Hamiltonian constraint that was generally used in LQC until recently (the so-called 
'/i -evolution'), the answer was in the negative. Certain quantum effects associated with that 
evolution could alter the classical predictions even in regimes in which the matter densities 
and space-time curvatures are completely tame. In the k=0 case, this severe drawback was 
overcome by the 'improved dynamics' of [3]. 

In the k=l models, the classical recollapse provides an excellent venue to test semi- 
classical viability because, as the analysis of section [VI] shows, both the classical recollapse 
and the quantum bounce are governed by the same condition: the vanishing of H 2 = (a/a) 2 . 
For the recollapse, an agreement with the classical theory requires that quantum geometry 
effects be negligible while the bounce can occur only if these effects dominate. At first sight 
then there appears to be a tension. In 0] , Green and Unruh analyzed this issue numerically 
in the same model as the one considered in this paper and concluded that the tension is real. 
More precisely, they used the then available Hamiltonian constraint of LQC 0] and found 
numerical evidence against the occurrence of recollapse. For large universes, the classical 
recollapse occurs when matter density and space-time curvatures are very small compared 
to the Planck scale. Therefore, a theory in which recollapse does not occur would contradict 
classical general relativity in a domain where there is every reason to expect its validity. 
Although Green and Unruh did not have access to a physical inner product or observables 
to arrive at a definitive interpretation of their results, they concluded that it is unlikely that 
one could find an interpretation in which such large deviations from the classical theory are 
appropriate. 

In this paper we have overcome important limitations pointed out by Green and Unruh, 
thereby completing the LQC program. We found that LQC does predict a recollapse and, 
furthermore, it occurs at the values of matter density and volume predicted by the classical 
theory. How did this strikingly different conclusion come about? We will conclude by first 
summarizing how the general criticisms of [6| were addressed and then discussing the issue 
of the recollapse. 

i) General Framework: Green and Unruh began by pointing out that much of the then 
successes of LQC arose from effective equations and the interpretation of quantum states 
had remained unclear because the physical inner product and observables had not been 
specified. They pointed out that, in particular, the issue of time had not been addressed 
explicitly and the implicit use of the scale factor as time has obvious problems in the k=l 
model. In this paper we showed that the use of the scalar field as emergent time is free 
of the difficulties associated with the multi-valued character of the scale factor in closed 
models. We constructed the physical sector of LQC in detail, including the physical inner 
product, Dirac observables and well-controlled semi-classical states. We then numerically 
solved the Hamiltonian constraint and calculated the expectation values and fluctuations 
of Dirac observables. In particular, we analyzed dynamics in full LQC, not in just in an 
effective approximation. 

ii) Recollapse: Our Hamiltonian constraint is quite different from that used in Thus, 
the very starting points of the two sets of numerical simulations are distinct. In particular, 
while Cgrav used in [6| was not self- adjoint, ours is. This has two consequences. First, since 
our operator 6 is self-adjoint on 7~t^ and we are guaranteed that it admits a complete set 
of eigenfunctions. Furthermore since its spectrum is discrete the eigenfunctions in the 
spectral family are normalizable. This immediately implies that they must decay for large \v\. 
Our numerical simulations showed that the decay is exponential for large which ensured 
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that physically appropriate wave packets would exhibit a recollapse. The numerical task of 
actually finding these eigenfunctions was delicate because the normalizable ones constitute 
a set of 'zero measure' among all eigenfunctions, i.e., because 'most' eigenfunctions diverge 
in at least one of the two asymptotic regimes v — > ±00. However, because is self-adjoint 
with appropriate properties, the existence of these eigenfunctions was ensured from the 
beginning. This was not the case for the analysis of p]. Indeed the principal argument there 
was that the authors found only exponentially growing eigenfunctions, signaling absence 
of recollapse. Secondly, for reasons discussed in SecJH it was important that we used the 
'improved dynamics' of [3] (rather than the older, L fi evolution'). Finally, the effective- 
equation-analysis shows how LQC manages to have both the quantum bounce and the 
classical recollapse. For, the modified Friedmann equation now has two roots: one a la 
classical general relativity at a low density causing the recollapse, and a new one near 
Planck density causing the bounce. The quantum geometry effects are small at the old, 
classical root but cause and thus dominate the new root. 
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APPENDIX A: INVARIANT FRAMES, EXPLICIT CHARTS AND HOLONOMY 

In this appendix we spell out conventions on the fiducial structures that are used in the 
main body of the paper and provide an explicit calculation of holonomy used in the definition 
of the Hamiltonian constraint. The discussion on conventions is somewhat detailed because 
there has been some confusion due to unfortunate typos in some of the standard literature 
on homogeneous cosmologies. 

1. Invariant frames 

In the k—1 case the underlying, spatial 3- manifold M is a 3-sphere S 3 . It is often conve- 
nient to identify it with the symmetry group SU(2) which acts on it simply and transitively. 
In what follows we will generally do so. Let us denote the symmetry vector fields on M 
by £f . All homogeneous isotropic tensor fields on M are invariant under diffeomorphisms 
generated by £f . In particular, these vector fields are the Killing vectors of all 3-metrics 
considered in this paper. 

Let us fix a basis t 1 in the Lie algebra su(2), satisfying rV J = |e y fcT fe — and denote 
by kij the metric on su(2) for which these t 1 constitute an orthonormal basis. In what follows 
the 'internal' or Lie-algebra indices will be lowered and raised using and its inverse k l K 
It follows in particular that := e mn kqmi(ljn is a 3-form on su(2) satisfying e^e*^ = 6. 

Recall that £77(2) admits a natural, left invariant, Lie-algebra valued, Cartan 1-form 
uj = g~ x dg. It naturally defines a co-frame u l a on M via: 



g ag = . uj = u) Ti 



(Al) 
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We will denote the dual frame by °e"; thus °e a i %o{ = 5\ and °e a i %j\ — 8%. From the definition 
of the natural left-invariant 1-form uo it follows that the co-frames 1o l a and the frames °e°: 
satisfy the relations: 

dV + \e\ jk °J A °u k and [%, °e 3 \ = e*/ °e k (A2) 

These will be our fiducial co-frames and frames on M. The 1-forms %j l a and the vector fields 
"ef are left invariant. Thus £f , the infinitesimal generators of left translations, Lie drag these 
fields: 

A* V = 0. ^ % = (A3) 

and satisfy the su(2) commutation relations = €ij k £k- The metric 

V:=°^°^^ (A4) 

on M will serve as the fiducial metric on M. By inspection £f are Killing fields of °q a b- Since 
they act simply and transitively on M, it follows that °q a b is of constant curvature. However, 
in contrast to one's initial expectations, it is the metric on a 3-sphere of radius a = 2 (rather 
than a = I). 12 

The volume of M with respect to °q a b is V = 2-ir 2 a 3 = 16tt 2 and its scalar curvature is 
°R = 6/a 2 = 3/2. It will be convenient for us to use the symbol V (rather than the numerical 
value 167r 2 ) to denote the fiducial volume of M. We will also set £ D = V^ 3 . In the standard 
charts used in textbooks, the components of °q ab can be expressed as: 



ds 2 Q = a 2 [dx 2 + sin 2 X (d6 2 + sin 2 9 dip 2 )] 
d '" +r 2 (d6 2 + sm 2 6d V 2 ) 



(A5) 



with a = 2, where x, ® £ (0,7r), <p G (0, 2n) and r = sin x G (0, 1). 

Finally, the fact that the choice a = 2 in the fiducial metric is geometrically natural can 
also be seen directly in terms of 3-sphere geometries without reference to SU(2) and the 
natural Cartan form thereon. Consider the Euclidean metric and 3-spheres §( a ) of radius a 
on M 4 . The natural action of the rotation group SO(4) on M 4 leaves each S( a ) invariant. It is 
the isometry group of the intrinsic metric on §( a ). These six Killing fields can be naturally 
divided into two SO (3) sub-Lie- algebras (resulting from self-dual and anti-self-dual 2- forms 
on IR 4 ). These are the right and left rotations and all three right rotations commute with all 
three left. In the natural chart x, y, z, w on M 4 , these six Killing fields Kf can be expressed 



This normalization is fixed by our choice that the natural left invariant 1-forms °ui l a be orthonormal. In 
the frame formalism it would be awkward and geometrically unnatural to work with rescaled °io % a . If one 
is interested only in metrics and not frames, on the other hand, one can just as easily work with an unit 
3-sphere. This is the usual choice in geometrodynamics. In the literature on cosmology the fiducial metric 
is sometimes chosen to have unit scalar curvature °R rather than unit radius; then the 3-sphere radius is 
a = l/v/6. 
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as: 

Kf = \ (xd y — yd x ± zdw =F wdz) 
K f = \ ( wd x - xd w ± zdy =f ydz) 

Kf — \ (ydw — wd y ± zdx =f xdz) (A6) 

so that [Kf,Kf\ = eij k K^ and [Kf,K~] = 0. Note that because of the first set of com- 
mutation relations, we do not have the freedom to rescale the Kf by a constant. The three 
vectors Kf in each set are mutually orthogonal. Now we can ask: On which S 3 is this basis 
orthonormal? The answer is: the 3-sphere with radius a = 2. On M, these are the basis 
and {#}. 



2. An explicit chart 

In the textbook treatments of SU(2), its elements are often written as 2 x 2 matrices: 

9 = ;f) where a = cos(/3/2) 

b = e-^ a ^ sin(/3/2) (A7) 

where < a < 2tt, < j3 < n and < 7 < Ait. (The ranges of a and 7 can be interchanged.) 
As is usual with angular coordinates, this chart breaks down at the poles a = 0, ( 5 = 0,7 = 
(which corresponds to the identity, I) and a = 0, (3 = ir, 7 = (which corresponds to —I). 
Nonetheless, as with the standard spherical polar coordinates, this chart is convenient for 
explicit calculations. 

In this chart, the left invariant co-frame °uj l a can be expressed as: 

°uj 1 = — cos 7d/3 — sin 7 sin (3da 
°u} 2 = sin 7d/3 — cos 7 sin [3da 

V 5 = d7 + cos/Ma (A8) 
and the left invariant frame e° as 

d sin 7 d cos (3 sin 7 <9 
ei = — COS7— : — -— — h 



d/3 sin /3 da sin /3 c?7 



„ d cos 7 9 cos (3 cos 7 9 
e 2 = sin 7— — h 



o. 



d{3 cos /9 da sin /3 ^7 
■e 3 = £. (A9) 



Finally the right invariant vector fields — which Lie drag our co-frame and the frame- are 
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given by: 

d cos (3 cos ad sin a d 

Si = C0SQ ^ ^5 — ^ _ + ^^^~ 

ap sin p aa sin p 07 

= — [cos a cos 7 — cos /3 sin a sin 7]%! + [cos a sin 7 + cos /3 sin a cos 7]°e 2 + [sin a sin /3]°e3 

. . d cos /3 cos a d cos a d 
4 2 = sin a— + 



dp sin /3 (9a sin (3 
= — [cos 7 sin a + cos (3 cos a sin 7]%! + [sin a sin 7 — cos /? cos a cos 7]°e 2 — [cos a sin /3]°e 3 

= sin/3sin7 °ei + sin/3 COS7 °e2 — cos/3°e3. (A10) 

Note that although the components of e" and diverge on the 2-torus f3 — 0, this is just 
an artifact of the failure of the chart there. All six vector fields are globally well-defined. In 
particular the expressions of as linear combinations of e° are globally well-behaved. 

From the expression of the 1-forms %f a we can readily compute the components of the 
fiducial metric: 

dS 2 = da 2 + d(3 2 + d 7 2 + 2dad 7 (All) 

The integral curves of the vector fields °e\ (as well as £3) are circles which provide a Hopf 
fibration of S 3 . The quotient M is topologically § 2 with the induced metric 

ds 2 = d(3 2 + sin 2 (3da 2 (A12) 

Thus (3, a are the standard polar coordinates on the quotient M and the induced metric on 
it is that of an unit 2-sphere. 



3. Holonomy 

The Hamiltonian constraint involves the field strength Fjf b of the SU(2)-connection 

K ■= c°ui (A13) 

In LQG, there is no operator corresponding to the connection itself; only the holonomy 
operators are well-defined. Therefore, we need to express the field strength as an appropriate 
limit of a suitable holonomy @, 0, Now in k=0 cosmologies, the right and left invariant 
vector fields on M coincide. Since they generate translations along the x, y, z directions, 
they commute. Hence to calculate the component F^ 2 := °e| °e\ F^ b of the field strength, we 
can construct a square □ by moving along z and y directions a distance p, (as measured 
by the fiducial metric °q a b), evaluate the holonomy along this loop, divide it by the area 
enclosed by the loop and take the limit as the loop shrinks. Because of homogeneity, F^ 2 
can be evaluated at any point of M and because of isotropy, F 3 fc 2 determines F^ b completely. 

In the present k—1 case we wish to follow a similar procedure. Again, F^ 2 can be evaluated 
at any point of M and fully determines F% b . However, now the vector fields °e\ and "e^ do 
not commute. Therefore, as explained in the main text, we can not form the desired loop 
□jj by moving along 4 segments of their integral curves. However, since the left invariant 
vector fields do commute with the right invariant vector fields £ 4 a , we can construct the 
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desired closed loop Qy using integral curves of °e" and 

For definiteness, let us start with the point (a = 0,/3 = tt/2, 7 = 0). At this point, 
£3 = and we can calculate the component °e\ °e b 2 F ab using the loop n 32 constructed 
from the composition of the following four segments: 

Segment 1: move from (0, tt/2,0) to (0, 7r/2, X£ q ) along the integral curve of 

°e 3 = d/d>y; 

Segment 2: move along (0, tt/2, X£ a ) to (X£ , tt/2, A£ ) along the integral curve of 
= 8 /da; 

Segment 3: move from (X£ Q , tt/2, \£ q ) to (X£ , tt/2, 0) along the integral curve of 
-°e 3 = -d/d>y; 

Segment 4: return from (X£ Q , tt/2,0) to the point of departure (0, tt/2,0) along 
the integral curve of £3 = —d/da. 

It is straightforward to calculate the holonomies along these four segments. One obtains 

h 1 = e XcT \ h 2 = e-MsinAn+cosAT^ /l3 = e -Acr 3 ^ = gAr 2 (A14) 

The holonomy ^□ 32 around the loop is just the composition /14/13/12/ii. A simple calculation 
yields: 

°es °4 *J = J™ Tr (^3 2 r fc ) = (c 2 - crf )*f (A15) 
Finally, using isotropy and homogeneity, one concludes: 

= - lim -L [sin 2 A(c - £ /2) - sin 2 (A£ /2)] e,/ °< ^ 

= C 2 (c 2 -c4)e,/K°^ (A16) 

The expression of the holonomy in the second step is used in the discussion of the Hamil- 
tonian constraint in section [IT] of the main text. 
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